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The  impact  on  a  prime  power  source  due  to  AC  resonant  charging  was 
investigated.  The  source  was  a  standard  DOD  15  KW,  400  Hz  alternator. 
Adverse  impact  on  the  alternator  was  almost  negligible  for  a  varying  load 
of  15  to  60  KW  and  a  frequency  change  of  380  to  475  Hz.  The  only  excep¬ 
tion  was  caused  by  the  pulser  diodes.  Their  turn-on  was  reflected  in 
the  alternator  voltage  waveforms.  An  indirect  approach  is  presented 
for  determining  resonance  and,  thus,  the  internal  reactance  contribution 
by  the  alternator.  This  value  was  found  to  be  slightly  less  than  that 
specified  for  the  leakage  reactance.  The  investigation  concludes  that 
resonance,  in  the  classical  steady  state  sense,  is  never  reached.  Even 
though  the  charging  circuit  has  a  transient  response,  the  alternator 
armature  reaction  is  approximately  fixed.  This  implies  that  the  alter¬ 
nator  is  in  a  near  steady  state  mode  of  operation. 


AC  RESONANT  CHARGING  FOR  INTERFACING 


PRIME  POWER  TO  PULSE  CONDITIONING  CIRCUITS 


I.  Introduction 


Background 

The  advent  of  directed  energy  weapons  has  spurred  the  development 
of  a  large  number  of  supporting  technologies.  Two  of  these  are  the 
power  source  and  power  conditioning.  The  development  of  these  techno¬ 
logy  areas  has  been  driven  by  several  factors.  Primary  among  these  are 
the  beam  weapon's  inefficiency  (requirement  for  high  power),  voltage  and 
current  waveshaping  (requirement  for  advanced  power  conditioning) ,  and 
dependable  interfacing  of  the  power  source  with  the  power  conditioning 
equipment  (requirement  for  lightweight,  low  volume  charging  of  the  inter¬ 
mediate  energy  store).  The  last  of  these  is  the  topic  of  concern  for 
this  research  effort. 

The  Army  and  Air  Force  are  both  concerned  with  the  interface  problem 
and  have  expended  considerable  resources  toward  its  resolution.  The 
Army's  need  stems  from  the  requirement  for  systems  that  are  ground  mobile. 
Thus,  all  equipment  in  this  general  category  must  be  relatively  light  and 
low  volume.  The  Air  Force  requirements  are  even  more  stringent  since  its 
systems  must  be  flyable. 

One  possible  way  to  resolve  the  interface  problem  is  by  the  use  of 
AC  resonant  charging.  This  method  employs  an  AC  voltage  source  to 
recharge  the  energy  storage  condenser  of  a  voltage-fed  network.  It  has 
been  the  subject  of  only  a  small  number  of  tests.  Thus,  there  is  only 
limited  information  (and  understanding)  concerning  the  interface  between 
the  alternator  and  the  charging  system  during  its  use. 


In  1976  Airesearch  Manufacturing  Company  submitted  their  findings 
for  the  AC  resonant  charging  tests  they  had  performed  on  a  modified 
alternator  (Ref  9).  One  of  the  secondary  objectives  of  the  program  was 
to  test  the  machine  and  its  loading  circuit  at  precise  resonant  condi¬ 
tions.  Although  this  objective  was  not  met,  the  machine  was  tested  at 
two  different  off-resonant  conditions.  This  study  was  sponsored  by  the 
Air  Force  Weapons  Laboratory  (AFWL) .  The  AFWL  also  sponsored  the  1978 
Power  Generation  Study  and  Test  Program  which  was  conducted  by  the  Mobi¬ 
lity  Equipment  Research  and  Development  Command  (MERADCOM) ,  Ft.  Belvoir, 
Virginia.  Among  other  conclusions,  it  found  that  including  all  of  the 
AC  resonant  charging  inductance  inside  the  alternator  was  not  only 
feasible  but  could  also  provide  an  overall  weight  saving  for  the  system. 
Most  of  the  charging  inductance  is  normally  included  in  a  high  leakage 
transformer. 

The  latest  effort  involving  AC  resonant  charging  is  a  program  being 
conducted  by  Clemson  University  and  sponsored  by  the  Aero  Propulsion  Lab¬ 
oratory  at  Wright-Patterson  AFB,  Ohio.  The  Clemson  investigators  are 
developing  computer  simulations  which  are  designed  to  predict  the  internal 
inductance  of  the  power  source  during  the  AC  charging  mode.  Thus,  the 
results  of  this  thesis  research  effort  should  serve  as  a  weighted  input 
to  either  refine  or  validate  the  Clemson  predictions. 

Problem  and  Scope 

The  major  problem  was  the  lack  of  a  full  understanding  of  the  inter¬ 
action  of  the  alternator  with  the  remainder  of  the  charging  circuit  during 
the  AC  resonant  charging  mode.  The  impact  on  the  alternator  needed  to  be 
determined  for  varying  loads  and  frequencies.  It  was  believed  that  this 
mode  of  operation  did  not  have  the  same  adverse  impact  on  the  alternator 
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as  did  DC  resonant  charging  (Ref  8:50).  However,  this  needed  to  be 
confirmed.  Also,  an  investigation  needed  to  be  made  to  determine  what 
actually  constitutes  resonance,  when  does  it  occur,  and  when  it  does 
occur  what  is  the  effective  internal  machine  reactance.  To  gain  a 
fuller  understanding  of  the  charging  system,  attention  also  had  to  be 
given  to  the  discharge  circuit,  system  grounding  (loops),  system  noise, 
pulse  triggering,  and  the  physical  layout  of  the  pulser.  Thus,  the  study 
was  primarily  concerned  with  the  charging  system  and  minimal  effort  was 
devoted  to  other  system  aspects  such  as  the  discharge  circuit. 

Approach 

The  first  step  in  this  investigation  was  a  review  of  AC  resonant 
charging  theory.  The  equatiors  derived  theoretically  were  then  consoli¬ 
dated  into  a  computer  program  called  EQUATE  that  developed  the  data  needed 
to  provide  plots  of  the  circuit  voltages  and  current.  Following  this  the 
Runge-Kutta  programs  RLCCHG  and  RLCDIS  were  written  to  integrate  the  dif¬ 
ferential  equations  of  the  charge  and  discharge  circuits,  respectively. 
These  latter  programs  solved  the  system  equations  for  all  frequencies 
while  program  EQUATE  was  applicable  only  at  the  circuit  resonant  fre¬ 
quency.  Once  these  theoretical  models  were  developed,  the  experimental 
effort  began.  The  system  was  configured  and  tested  at  15,  30,  45,  and 
briefly,  at  60  kilowatts.  Special  tests  were  also  made  in  order  to 
determine  resonance  and,  thus,  the  effective  machine  reactance. 

Sequence  of  Presentation 

A  detailed  description  of  AC  resonant  charging  theory  is  presented 
in  Chapter  II.  Chapter  III  describes  the  experimental  results  and  their 
analyses.  The  project  conclusions  and  recommendations  for  future  study 
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are  presented  in  the  final  chapter. 

There  are  also  five  appendices.  Appendix  A  presents  a  detailed 
development  of  the  half  cycle  charging  equations.  Appendix  B  provides 
a  listing  of  program  EQUATE.  Appendix  C  includes  the  operating  instruc¬ 
tions  for  developing  programs  which  will  use  the  subroutine  RKFOUR.  It 
also  provides  a  listing  of  RKFOUR,  RLCCHG,  and  RLCDIS.  Appendix  D  con¬ 
tains  a  list  of  the  test  equipment.  The  last  appendix  gives  the  results 
of  the  leakage  transformer  short-circuit  test. 
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II.  Theory  of  AC  Resonant  Chargin; 


While  the  theory  of  AC  resonant  charging  has  been  known  for  at 
least  the  last  four  decades  it  has  yet  to  receive  thorough  publication 
coverage.  The  "grandfather  text"  is  Pulse  Generators  (Ref  2)  which  was 
edited  by  Glasoe  and  Lebacqz  shortly  after  World  War  II.  The  following 
material  is  a  brief  synopsis  taken  from  this  text. 

The  Pulser 

The  circuit  that  delivers  the  pulses  of  energy  to  a  load  is  referred 
to  as  a  "line-type  pulser."  This  type  of  pulser  discharges  all  of  its 
energy  during  each  pulse.  The  scope  of  this  presentation  does  not  cover 
the  "hard-tube  pulser"  which  delivers  only  a  small  fraction  of  its  stored 
energy  into  the  load  during  a  single  pulse.  Generally,  the  discussion  of 
a  pulser  can  logically  be  divided  into  two  categories,  the  charge  and  dis¬ 
charge  circuits.  This  is  possible  since  the  charging  of  the  energy-storage 
component  of  the  pulser  requires  a  time  that  greatly  exceeds  the  discharge 
time.  The  discharge  circuit  will  be  introduced.  However,  this  presenta¬ 
tion  will  be  directed  primarily  at  the  charging  circuit. 

The  Discharge  Circuit 

The  basic  discharge  circuit  of  a  pulser  can  be  represented  schemati¬ 
cally  as  shown  in  Figure  1.  If  the  energy  for  the  pulse  is  stored  in  the 
electrostatic  field  of  a  capacitor (s) ,  then  the  unit  is  a  "voltage-fed 
network"  as  opposed  to  a  "current-fed  network"  which  stores  its  energy 
in  a  magnetic  field.  This  energy  storage  device  is  often  referred  to  as 
a  "pulse-forming  network"  or  PFN.  "When  the  load  impedance  is  equal  to 
the  characteristic  impedance  of  the  network,  assuming  the  switch  to  have 
negligible  resistance,  all  the  energy  stored  in  the  network  is  transferred 
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Figure  1.  Tvpical  Discharge  Circuit  of  a  Pulser 
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to  the  load,  leaving  the  condensers  in  the  network  completely  discharged" 
(Ref  2:8). 

A  simplified  schematic  for  the  pulser  discharge  circuit  built  for 
the  MERADCOM  facility  is  shown  in  Figure  2.  A  single  capacitor  is  used 
as  the  energy  storage  device  and,  thus,  the  circuit  does  not  truly  con¬ 
tain  a  PFN.  The  load  is  an  eleven  ohm  resistor  and  a  thyratron  is  used 
for  switching.  A  small  inductance  has  been  placed  in  the  series  loop  in 
order  to  drive  the  thyratron  anode  voltage  to  a  slightly  negative  value 
which  assures  turn-off  of  the  switch.  Therefore,  all  of  the  available 
energy  is  not  discharged  to  the  load.  A  very  small  amount  is  returned 
and  negatively  charges  the  energy  storage  capacitor,  This  process,  how¬ 
ever,  does  not  consume  energy  but  rather  slightly  reduces  the  amount 
which  can  be  used  by  the  load. 

The  Charging  Circuit 

Resonant  charging  of  a  voltage-fed  network  in  a  line-type  pulser  can 
be  accomplished  in  several  ways.  For  example,  one  of  the  best  known 
methods  is  to  use  a  DC  power  supply.  This  approach  was  studied  by  Lt. 
Jaime  Silva  who  noted  that  it  presented  a  number  of  problems,  namely, 
"vibration  of  the  generator  at  the  pulse  repetition  rate,  irregular 
voltage  and  current  waveforms  in  the  generator,  and  lower  performance 
of  the  DC  resonant  circuit  due  to  input  voltage  sag"  (Ref  8:50).  How¬ 
ever,  a  system  based  upon  this  method  is  relatively  simple  in  its  con¬ 
struction  and  triggering  features. 

Resonant  charging  may  also  be  realized  using  an  AC  power  supply. 
Though  more  complicated,  this  approach  has  a  number  of  advantages  over 
DC  resonant  charging.  As  concluded  by  Lt.  Silva,  this  latter  charging 
scheiRT  is  not  expected  to  cause  a  negative  impact  upon  the  power  source. 
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However,  it  does  require  that  the  AC  source  frequency  be  integrally 
related  to  the  pulse  recurrence  frequency  (PRF) .  This  is  necessary 
since  the  charging  period  must  terminate  at  a  current  zero  crossing. 

As  there  are  only  two  zero  crossings  per  cycle,  the  PRF  must  not  be 
greater  than  twice  the  AC  frequency.  The  relationship  between  the  PRF 
and  the  AC  frequency  is  given  by  the  following  equation. 

?RF  -  |  fAC  (1) 

where 

K  =  1,2... 

f^c  =  impressed  AC  frequency 

A  value  of  K  =  1  corresponds  to  half  cycle  charging  and  K.  =  2  to  full 
cycle  charging.  If  the  PRF  is  twice  the  AC  frequency  then  the  network 
voltage  attains  a  value  approximately  tt/2  times  the  peak  AC  voltage. 
Though  this  fixed  frequency  requirement  allows  little  flexibility  in  the 
altemator-pulser  interface,  the  deficiency  is  offset  somewhat  by  the 
advantage  of  being  able  to  control  the  pulse  power  output  by  varying  the 
alternator  field  current.  Another  advantage  is  that  this  system  permits 
a  net  saving  in  weight  and  size  relative  to  the  DC  system.  Its  block 
diagram  is  shown  in  Figure  3. 

There  are  a  number  of  requirements  levied  upon  this  circuit.  Among 
them  is  the  need  for  high  efficiency.  This  normally  eliminates  the  use 
of  a  resistor  for  the  charging  element  since  its  inherent  efficiency  can 
never  exceed  50  per  cent.  The  charging  element  must  allow  the  discharge 
of  the  network  at  a  peak  voltage.  It  must  also  be  capable  of  isolating 
the  source  from  the  time  the  switch  closes  until  it  deionizes  (for  a 
gaseous  switch) .  For  example,  if  the  network  is  connected  directly  to 


Figure  3.  Block  Diagram  for  the  Complete  Pulser  Circuit 


Figure  4.  Schematic  Diagram  of  the  MERADCOM  Pulser 
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the  high  voltage  terminals  of  a  transformer  and  discharged  at  a  voltage 
peak,  the  network  begins  to  recharge  immediately,  the  switch  fails  to 
re-open,  and  the  transformer  is  short-circuited  for  the  remainder  of 
the  half  cycle. 

This  problem  may  be  avoided  by  placing  a  hold-off  diode  between  the 
high  voltage  transformer  and  the  energy  storage  network.  Figure  8  shows 
such  an  arrangement.  The  network  (capacitor)  may  now  be  discharged  during 
the  following  half  cycle  which  has  a  polarity  opposite  to  that  of  the 
charging  half  cycle.  This  provides  stable  circuit  action  and  eliminates 
the  problem  associated  with  the  continual  build-up  of  the  charging  voltage 
and  current  should  the  discharge  circuit  switch  fail  to  operate  (close) . 

As  shown  in  Figure  8,  all  phases  have  two  charging  circuits  and  each  of 
these  contains  a  hold-off  diode.  One  circuit  network  is  charged  during 
a  positive  half  cycle  and  the  other  is  charged  during  a  negative  half 
cycle.  The  transformer  utilization  factor  is  greatly  improved  by  this 
approach. 

The  MERADCOM  Pulser 

These  advantages  are  incorporated  in  the  circuit  of  Figure  4.  This 
circuit  is  the  basis  for  the  MERADCOM  pulser.  It  provides  a  more  detailed 
account  of  one  of  the  phase  networks  of  Figure  8.  A  description  of  its 
operation  will  now  be  presented. 

During  the  positive  half  cycle  capacitor  Cl,  the  energy  storage 
device,  is  charged.  The  voltage  on  Cl,  V^,  reaches  its  peak  when  the 
source  (alternator)  voltage  and  charging  current  pass  through  zero.  VC1 
remains  at  this  peak  value  for  a  brief  time  into  the  negative  half  cycle. 
When  the  source  voltage  reaches  a  specified  negative  value  a  trigger 
pulse  is  sent  to  the  thyratron  in  the  Cl  discharge  circuit.  The  thyratron 


fires  and  discharges  the  stored  energy  into  the  11  ohm  resistive  load. 

As  the  discharge  current  tends  to  zero, the  voltage  on  the  small  inductor 
at  the  switch  anode  reverses  polarity.  This  forces  a  slight  negative 
charge  on  Cl  which  assures  proper  turn-off  of  the  thyratron.  During  the 
negative  half  cycle  the  source  charges  C4  and  the  same  chain  of  events 
occur  as  described  for  Cl. 

General  Approach  to  the  Theoretical  Analysis 

Independent  theoretical  analyses  were  performed  on  both  the  discharge 
circuit  and  the  charge  circuity  of  the  pulser.  A  number  of  computer  pro¬ 
grams  were  developed.  These  divide  into  twogene?Si--eategories  as  implied 
in  the  introduction.  The  first  of  these  was  simply  an  application  of  the 
AC  resonant  charging  equations  derived  in  Appendix  A.  These  equations 
were  used  to  show  how  the  circuit  current  and  voltages  varied  over  a 
half  cycle  of  charging.  However,  their  use  was  applicable  only  to  a  cir¬ 
cuit  operating  at  the  resonant  frequency.  This  program  of  equations  is 
listed  in  Appendix  B  as  Program  EQUATE. 

The  second  program  category  was  based  upon  the  use  of  a  Runge-Kutta 
integration  subroutine  called  RKFOUR  which  was  developed  by  B.  D.  Weathers. 
Programs  were  developed,  following  his  suggested  format,  that  would  inte¬ 
grate  the  differential  equations  for  a  number  of  cases.  The  most  useful 
of  these  were  Program  RLCCHG  and  Program  RLCDIS  which  integrated  the 
charge  and  discharge  circuit  equations,  respectively.  Program  RLCDIS 
was  obtained  by  slightly  modifying  Program  RLCCHG.  Both  of  these  programs, 
as  well  as  RKFOUR  and  instructions  for  its  use,  are  listed  in  Appendix  C. 
They  each  integrate  their  respective  circuit  equations  over  a  half  cycle. 
Program  RLCCHG  was  found  to  be  particularly  useful  since  it,  in  contrast 
to  Program  EQUATE,  is  capable  of  integrating  the  charging  circuit  equations 
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As  noted  earlier,  the  discharge  circuit  for  the  MERADCOM  pulser  can 


be  represented  by  the  schematic  shown  in  Figure  2.  The  differential 
equations  for  this  circuit  were  integrated  using  Program  RLCD1S  with  the 
resistance  of  the  switch  neglected.  The  current  and  voltage  curves  shown 
in  Figures  5  and  6  resulted  from  the  selection  of  a  0.43  yf  capacitor, 
an  11  ohm  load  resistor,  and  a  29  yh  inductor.  The  capacitor  was 


TIME 


Figure  5.  Pulser  Discharge  Current  Using  Program  RLCDIS 
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Figure  6.  Discharge  Circuit  Capacitive  and  Inductive  Voltages 
Using  Program  RLCDIS 
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initially  charged  to  -6000  volts.  This  value  was  read  from  the  oscillo¬ 
grams  for  a  load  of  approximately  20  KW.  As  seen  in  the  curves,  at  the 
time  of  the  current  zero  crossing  the  capacitor  voltage  has  reversed  its 
polarity.  This  reverse  charge  forces  the  thyratron  switch  to  turn  off. 
Thus,  the  discharge  circuit  is  isolated  from  the  charging  circuit  through¬ 
out  the  interpulse  period. 

Theoretical  Analysis  of  the  Charging  Circuit 

Figure  7  shows  the  simplified  schematic  for  the  circuit  used  in  half 
cycle  AC  resonant  charging.  The  more  detailed  three  phase  representation 
of  the  MERADCOM  circuit  is  shown  in  Figure  8. 


Figure  7.  Simplified  Circuit  for  Half  Cycle  Charging 


The  mathematical  analysis  for  the  circuit  of  Figure  7  is  given  in 
Appendix  A.  The  charging  inductance,  Lc,  is  made  up  of  the  internal 
machine  inductance,  L  ,  and  the  large  transformer  leakage  inductance,  L  . 
The  charging  resistance,  Rc,  is  the  series  combination  of  the  source 
internal  resistance,  the  transformer  resistance,  and  the  load  resistance. 
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is  the  fixed  capacitance  of  the  energy  storage  device. 

From  Appendix  A  it  was  determined  that  Cjj  charged  up  to  its  maximum 
voltage  when  the  power  source  phase  angle  was  equal  to  -tt/2.  Using  this 
value  and  recognizing  that  q^(0)  =  C^vN (0)  the  following  expressions  for 
the  charging  current  and  the  circuit  voltages  may  be  obtained  (see 
Appendix  A) 

£0  t 

a 

ic(t)  =  QCNEaa)a((l  -  e"  2Q)sinoJat 


6J  t 
a 


v  (0) _ 

~QE  e  2(3  sinuat} 


oj  t 

a 


vN(t)  =  Q{  (e  2Q  -  l)costoat 


_  VN(0>  i  i  v  (0) 

2Q  -  cos^t  +  (^  +  ^ 


(2) 


(3) 
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where 

.  * 

<o  L 

Q  •  — ^  c  ,  the  quality  factor  of  the  circuit 
c 

E&  =  peak  emf  generated  inside  the  alternator  (volts) 

=  the  natural  resonant  frequency  (radians/sc - ond) 
vN(0)  =  the  initial  capacitor  voltage  (volts) 
v^(t)  =  the  voltage  across  the  inductor  (volts) 

Equations  ( 2) — ( 4)  may  now  be  plotted  versus  time.  When  this  is  done 
with  Program  EQUATE  the  curves  in  Figures  9  and  10  result.  Figure  9 
closely  resembles  the  typical  tsint  curve.  This  form  is  evident  in  the 
simplified  expression  for  the  charging  current.  For  example,  with 
1^(0)  =  0,  $  =  — tt /2 ,  and  qN(0)  •  CNvfJ(0),  equation  (A-22)  becomes 
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ohms,  C„  =  147  pf,  v.,(0)=-5.4  volts,  E  =  150  volts,  and  L  =  0.858  mh. 

N  N  ’a  c 

These  values  were  selected  early  in  this  study.  Their  substitution  into 
equation  (5)  yields 

ic(t)  =  (87413t  +  1.99)sin2802t  (6) 

The  Results  section  will  show  that  the  above  values  are  not  very  different 
from  the  final  values  which  were  either  measured  or  calculated. 

With  the  exception  of  the  first  22  microseconds,  the  complete  half 
charging  cycle  (1.12  msec)  is  dominated  by  the  first  term  in  equation  (6), 
that  is,  the  O^tsinw^t  term.  Since  the  charging  voltage  v^(t)  is  a  func¬ 
tion  of  the  integral  of  the  current  it  will  reach  its  maximum  value  at 
the  end  of  the  charging  half  cycle  (it  radians) .  Again,  the  purpose  of 
the  diode  in  Figure  7  is  to  block  the  circuit's  attempt  to  lower  this 
voltage  after  it  has  reached  its  peak  value.  Without  the  diode  the 
current  would  go  negative  at  the  end  of  the  half  cycle.  This  would  cause 
a  reverse  polarity  on  the  capacitor  which,  at  the  end  of  a  full  cycle, 
would  be  twice  the  half  cycle  value.  This  is  known  as  full  cycle  charging. 
Without  proper  controls  the  voltage  could,  as  explained  earlier,  build  up 
to  a  dangerous  level,  especially  should  the  discharge  switch  fail  to 
operate.  With  half  cycle  charging  this  disadvantage  is  eliminated. 
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III.  Results  and  Discussion 


The  circuit  of  Figure  8,  with  the  discharge  switch  branch  included, 
was  set  up  and  testing  was  begun.  However,  the  circuit  was  later  par¬ 
tially  disconnected  to  confirm  the  resistance,  inductance,  and  capaci¬ 
tance  values  of  the  major  components.  These  components  are  listed  in 
some  detail  in  Appendix  D.  The  per  phase  DC  resistance  of  the  alterna¬ 
tor  was  found  to  be  approximately  0.18  ohms.  This  value  was  multiplied 
by  a  factor  of  1.5  to  account  for  the  resistive  change  due  to  loss 
mechanisms  such  as  hysteresis  and  eddy  currents  (Ref  4:255;  5:227-237). 

A  short  circuit  test  of  the  C  phase  high  leakage  transformer  was  then 
conducted.  The  results  of  this  test  are  recorded  in  Appendix  E  where 
all  values  are  referred  to  the  low  voltage  side  of  the  transformer.  The 
resistance  and  inductance  (at  400  Hz)  were  found  to  be  0.1427  ohms  and 
0.690  mh,  respectively.  Also,  a  capacitance  of  0.415  pf  was  measured 
for  one  of  the  energy  storage  (charging)  capacitors.  Since  these  tests 
were  not  performed  initially,  some  of  the  earlier  theoretical  analyses 
(plots,  for  example)  were  based  on  slightly  different  nominal  component 
values . 

System  Grounding  and  Noise 

During  the  early  stages  of  testing  the  pulser  system  was  found  to 
be  inherently  noisy.  It  was  this  problem  that  restricted  Lt  Silva's 
operation  of  the  system  to  a  low  power  level.  His  report  indicates  that 
multiple  trigger  signals  prevented  the  build-up  of  voltage  on  the  charging 
capacitors  and  that  this  malfunction  was  caused  by  the  electronic  logic 
circuit  (Ref  8) .  He  did  not  have  sufficient  time  to  diagnose  this  prob¬ 
lem.  However,  after  his  departure  the  MERADC0M  engineers  discovered  that 
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one  of  the  wires  which  fed  the  trigger  pulse  to  a  thyratron  switch  was 
positioned  too  close  to  one  of  the  discharge  circuit  cables.  Such  close 
proximity  allowed  each  discharge  in  the  cable  to  cause  a  corresponding 
noise  pulse  in  the  thyratron  trigger  circuit.  Thus,  the  noise  from  one 
discharge  circuit  caused  the  inadvertent  discharge  of  another  circuit. 
The  charging  circuit  inductance  isolated  the  power  source  during  this 
positive  half  cycle  transient.  Such  isolation  is  normally  needed  only 
during  the  negative  half  cycle  since  the  storage  capacitors  are  dis¬ 
charged  during  this  period  of  time.  The  hold-off  diodes  perform  this 
isolation  function. 

The  physical  layout  of  the  pulser  complemented  the  noise  problem  by 
causing  a  number  of  ground  loops.  Initially,  there  were  ground  ties  at 
the  front,  rear,  and  sides  of  the  pulser  structure.  With  certain  mea¬ 
surements  this  caused  no  problems.  However,  some  measurements  could  not 
be  taken  simultaneously  with  the  dual  trace  scope.  For  example,  it  was 
not  possible  to  initially  measure  the  charging  current  and  charging 
voltage  at  the  same  time.  The  charging  current  measurement  was  obtained 
by  use  of  a  current  transformer  which  was  grounded  at  the  right  end  of 
the  structure.  The  ground  for  the  charging  voltage  was  located  at  the 
front  of  the  structure.  When  these  two  grounds  were  tied  together  at 
the  scope  inputs,  the  display  of  the  charging  current  waveform  was  very 
distorted.  The  system  was  improved  somewhat  by  either  removing  (as  in 
the  case  of  the  charging  current  ground)  or  relocating  (as  with  the 
charging  voltage  ground)  certain  ground  connections. 

It  quickly  became  evident  that  considerable  improvement  could  have 
been  made  in  the  shielding,  grounding,  and  physical  layout  of  the  pulser. 
It  would  be  difficult  to  correct  these  deficiencies  now  that  the  pulser 


has  been  constructed.  Hindsight  suggests,  as  usual,  that  they  should 
have  been  made  during  the  design. 

The  Discharge  Circuit 

Though  the  discharge  circuit  was  not  studied  thoroughly,  its  cur¬ 
rent  waveform  was  viewed  on  the  oscilloscope  in  order  to  make  a  theore¬ 
tical  comparison.  By  integrating  the  equations  for  the  simplified  dis¬ 
charge  circuit  of  Figure  2  the  discharge  current  waveform  of  Figure  5 
was  obtained.  The  small  thyratron  turn-off  inductance  of  the  circuit  was 
initially  measured  to  be  approximately  51  Uh.  This  caused  the  post-dis¬ 
charge  voltage  on  the  capacitor  to  reach  an  abnormally  high  level  (approx¬ 
imately  800  volts) .  Reducing  the  inductance  to  29  Uh  by  tap  adjustment 
caused  this  voltage  to  decrease  to  a  value  less  than  300  volts.  Since 
this  value  of  voltage  was  quite  sufficient  for  turning  the  thyratron  off, 
it  was  used  as  an  integral  coefficient  in  the  development  of  the  curve  in 
Figure  5.  This  curve  may  be  compared  with  the  measured  discharge  current 
curve  of  Figure  11. 

These  waveforms  are  quite  similar  and  some  of  their  differences  can 
be  readily  explained.  For  example,  the  calculated  current  peak  approaches 
340  amps  whereas  the  measured  current  peak  is  only  300  amps.  Some  reasons 
for  this  difference  may  be  the  use  of  a  0.43  uf  capacitor  in  the  calcula¬ 
tions  while  its  measured  value  was  found  to  be  0.415  Uf,  the  neglecting 
of  all  circuit  resistance  except  the  load,  and  the  inclusion  of  only  the 
thyratron  turn-off  inductance.  Each  of  these  modeling  assumptions  tend 
to  increase  the  calculated  current  peak.  Also,  the  computer  program 
selection  of  an  initial  capacitor  voltage  of  6000  volts  was  based  upon 
a  scope  reading  that  could  well  be  off  by  5%. 
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Figure  11.  Pulser  Discharge  Current  (50  amps/div,  R  =  11  ohms, 
C  =  0.415  yf,  L  =  29  yh,  =  6000  volts) 


Figure  12.  Primary  Current  and  Voltage  at  the  Alternator  Terminals 
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The  second  difference  in  the  waveforms  is  in  the  time  to  the  cur¬ 


rent  zero  crossing.  The  calculated  waveform  shows  a  half  period  of  only 
15  ysec  while  the  experimental  value  approaches  19  ysec.  Correcting  the 
capacitance  to  the  smaller  value  tends  to  decrease  the  calculated  half 
period  even  further.  However,  an  increase  in  the  circuit  resistance  and 
inductance  due  to  the  various  components,  such  as  the  thyratron,  tend  to 
increase  the  half  period.  Thus,  a  correction  in  the  circuit  resistance 
and  inductance  appears  to  override  the  slight  change  in  capacitance  in 
regard  to  the  discharge  waveform. 

The  Charging  Circuit 

Alternator  Terminal  Measurements  of  Voltage  and  Current.  Various 
measurements  were  made  at  both  primary  and  secondary  points  for  the  cir¬ 
cuit  of  Figure  8.  The  current  and  voltage  waveforms  at  the  alternator 
terminals  were  needed  to  help  determine  the  impact  of  resonant  loading 
on  the  power  source.  The  oscillogram  of  Figure  12  shows  the  primary 
current  and  voltage  at  the  alternator  terminals  with  the  system  operating 
at  a  speed  of  2250  RPM  (450  Hz)  and  excited  to  an  output  power  level  of 
45  KW.  All  system  measurements  were  taken  with  a  constant  DC  voltage 
applied  to  the  exciter  field  since  a  satisfactory  voltage  regulator  was 
not  available  for  this  function. 

Though  some  differences  may  be  noted,  the  general  appearance  of  the 
waveforms  is  sinusoidal.  This  was,  in  particular,  expected  for  the 
voltage  waveform.  It  contains  a  number  of  abrupt  changes  which  are  due 
to  the  turning  on  of  the  diodes  in  the  pulser.  Note  that  there  are  six 
transients  per  cycle.  This  number  corresponds  to  the  number  of  diodes. 
The  two  largest  transients  occur  near  the  C  phase  voltage  zero  crossing 
and  reflect  the  turn  on  of  the  two  diodes  in  this  particular  phase.  Any 
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time  a  diode  in  one  phase  turns  on  it  disturbs  the  resultant  armature 
reaction  flux  and,  thus,  the  resultant  air  gap  flux.  This  disturbance 
affects  all  phases,  especially  the  phase  containing  the  diode  that  turned 
on.  One  other  point  to  be  noted  is  that  the  voltage  does  not  experience 
any  sagging  as  the  current  rises  to  a  peak.  This  was  one  of  the  problems 
reported  by  Lt  Silva  in  his  investigation  of  DC  resonant  charging  (Ref  8: 
32-36) . 

The  current  waveform  in  Figure  12  presents  the  typical  tsint  shape. 
Recall  from  equation  (4)  that  the  theory  predicted  this  type  of  waveform. 
Though  not  pictured,  the  waveforms  for  15,  30,  and  60  KW  had  the  same 
general  appearance  as  those  shown  in  Figure  12. 

Pulser  Measurements.  Figure  13  shows  the  charging  current  and  volt¬ 
age  waveforms  for  one  of  the  two  charging  circuits  of  Figure  4.  The  mea¬ 
surements  were  made  at  a  point  between  the  diode  and  the  0-5  amp  ammeter. 
The  machine  was  delivering  approximately  20  KW  at  2200  RPM  (440  Hz) .  As 
in  Figure  12,  the  half  cycle  charging  current  waveform  displays  the  basic 
tsint  shape.  The  general  form  of  this  curve  compares  favorably  with  the 
calculated  charging  current  curve  in  Figure  9.  However,  these  curves 
cannot  be  compared  on  a  magnitude  basis  since  a  value  for  the  driving 
voltage  is  assumed  in  Figure  9  but  is  unknown  (and  cannot  very  easily  be 
determined)  for  the  oscillogram. 

The  charging  voltage  waveshape  in  the  oscillogram  closely  resembles 
the  VC  curve  in  Figure  10.  However,  the  calculated  voltage  values  in 
Figure  10  have  been  referred  to  the  transformer  primary.  This,  plus  the 
problem  concerning  the  unknown  internal  emf,  means  that  the  two  voltage 
waveforms  cannot  be  compared  on  a  quantitative  basis. 
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Figure  13.  Charging  Current  and  Voltage 


Figure  14.  Initial  Rise  of  the  Charging  Current 


The  charging  voltage  and  current  waveshapes  for  various  frequencies 
and  power  levels  present  a  number  of  unexpected  deviations.  These  anom¬ 
alies  will  now  be  discussed. 

Charging  Current  Anomalies.  One  of  the  irregularities  associated 
with  the  charging  current  waveform  occurs  at  the  end  of  the  half  cycle. 
Figure  13  shows  that  the  current  passes  through  zero  and  briefly  stays 
negative.  A  close  examination  of  the  charging  voltage  reveals  that  it 
decreases  slightly  at  this  same  time.  This  implies  that  the  current 
reverses  direction  and  reduces  the  charge  on  the  capacitor  by  a  minute 
amount.  These  actions  are  possible  only  because  the  hold-off  diode  is 
a  combination  of  resistors,  capacitors,  and  diodes;  that  is,  it  is  a 
stacked  rectifier.  Its  components  are  listed  in  Appendix  D.  While  the 
charging  current  is  positive  the  compensating  capacitors  in  the  rectifier 
are  fully  charged  in  one  direction.  As  the  source  (high  voltage  trans¬ 
former  terminals)  begins  to  change  its  polarity  the  large  voltage  on  the 
charging  capacitor  forces  an  immediate  voltage  reversal  on  the  compen¬ 
sating  capacitors.  This  also  causes  the  circuit  current  to  briefly  go 
negative.  This  current  removes  charge  from  the  storage  capacitor  to  the 
compensating  capacitors  which,  when  fully  charged,  reduce  the  current 
back  to  zero.  The  amount  of  charge  involved  in  this  transfer  is  so 
small  that  it  is  barely  detectable  on  the  voltage  oscillogram. 

The  charging  current  waveform  also  has  an  irregularity  at  the 
beginning  of  the  half  cycle.  The  initial  rise  in  the  charging  current 
of  Figure  13  is  shown  in  Figure  14.  Again,  the  hold-off  diode  is 
believed  to  be  responsible  for  this  deviation.  The  vertical  transient 
of  approximately  0.3  amps  is  difficult  to  explain  since  the  charging 
circuit  contains  such  a  large  charging  inductance.  This  inductance  of 


approximately  0.3  henries  (referred  to  the  secondary)  should  eliminate 
any  high  di/dt  component.  The  rounded  portion  of  the  initial  current 
rise  seems  to  be  caused  by  the  resonating  effect  of  the  charging  induct¬ 
ance  with  the  compensating  capacitors  in  the  rectifier  stack.  These 
components  produce  a  response  which  reaches  a  peak  in  approximately  35.7 
ysec.  Though  this  value  is  quite  different  from  the  50  ysec  oscillogram 
time,  it  seems  to  be  a  reasonable  method  for  analyzing  the  current 
initiation. 

Charging  Voltage  Anomaly.  The  charging  voltage  curve  in  Figure  13 
looks  almost  like  the  ideal  textbook  curve.  However,  the  initial  tests 
resulted  in  the  curve  shown  in  Figure  15.  The  charging  voltage  begins 
at  zero  or  some  negative  value.  After  a  brief  rise  it  suddenly  experi¬ 
ences  a  transient.  A  blowup  of  a  similar  transient  observed  in  a  later 
test  is  shown  in  Figure  16. 

There  are  two  major  reasons  for  the  difference  in  the  voltage  curves 
of  Figures  13  and  15.  The  alternator  speed  and  thyratron  hold-off  induct¬ 
ance  for  the  transient-plagued  curve  of  Figure  15  were  2000  RPM  and  51  yh, 
respectively.  The  smooth  curve  of  Figure  13  was  obtained  at  a  speed  of 
2200  RPM  and  the  inductance  was  reduced  to  29  yh.  The  reduction  in  the 
inductance  was  primarily  made  in  order  to  lower  the  reverse  voltage  on 
the  charging  capacitor  at  the  end  of  each  discharge.  A  reduction  in  the 
charging  voltage  transient  turned  out  to  be  a  positive  side  effect  of 
this  change.  Though  not  as  evident,  the  transient  reappeared  at  higher 
power  levels.  Thus,  the  changes  in  frequency  and  inductance  helped  but 
did  not  eliminate  the  source  of  the  problem. 

Since  the  magnitude  of  the  transient  depended  upon  the  value  of  the 
thyratron  turn-off  inductance,  the  possibility  of  an  interaction  between 
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Figure  15.  Initial  Charging  Voltage  Curve  with  Transient 


Figure  16.  Blowup  of  Charging  Voltage  Curve  Transient 
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the  charge  and  discharge  circuits  of  the  phase  was  investigated.  The 
charging  voltage  curve  for  one  circuit  of  the  phase  was  observed  after 
disconnecting  the  external  trigger  input  to  the  thyratron  grid  in  the 
second  circuit.  This  effectively  removed  the  second  circuit  from  the 
phase.  When  this  was  done  the  transient  was  eliminated. 

Thus,  it  appears  that  the  problem  stems  from  voltage  spikes  produced 
by  the  thyratron  turn-off  inductance  in  each  discharge  circuit.  In  order 
to  analyze  the  interaction  between  the  circuits  the  schematic  diagram  of 
the  pulsor  has  been  simplified  as  shown  in  Figure  17.  The  circuit  analysis 
may  begin  by  assuming  that  Cl  has  been  fully  charged  and  is  awaiting  dis¬ 
charge  through  thyratron  Tl.  The  charging  of  C4  has  just  begun.  Thus, 

CRl  is  turned  off  and  CP.4  is  turned  on.  If  Tl  is  now  triggered  to  close. 


Figure  17.  Simplified  Pulser  Schematic  for  Studying 
the  Charge/Discharge  Interaction 


the  discharge  current  through  it  causes  the  voltage  on  L^(VL^)  t0  S° 
positive  and  then  negative.  As  the  current  approaches  zero  this 
reversed  polarity  voltage  reaches  a  magnitude  that  exceeds  the  slow 
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rising  voltage  on  capacitor  C4.  (It  also  exceeds  the  source  voltage. 
However,  the  charging  inductance,  LC,  completely  isolates  the  source 
from  this  transient.)  When  the  voltage  on  C4  is  exceeded,  CR1  turns  on. 
The  charge  which  would  have  gone  into  building  up  additional  reverse 
bias  on  Cl  is  now  shunted  onto  C4.  The  current  path  for  this  charge  is 
LI,  Tl,  R4,  C4,  CR4,  and  CR1.  When  the  discharge  current  passes  through 
zero,  switch  Tl  opens,  V  ^  drops  back  to  zero,  and  the  transient  ends. 

From  the  transient  standpoint  it  would  be  sufficient  to  remove  LI 
and  L4  from  the  circuit.  However,  they  are  required  in  order  to  assure 
the  turn-off  of  the  thyratrons.  Thus,  for  this  particular  pulser  the 
symptom,  rather  than  the  cause,  must  be  addressed.  One  approach  is  to 
incorporate  an  SCR  in  series  with  each  rectifier  stack.  The  trigger 
timing  should  allow  the  SCR  to  turn  on  only  during  the  charging  half 
cycle.  Therefore,  even  if  the  transient  should  forward  bias  the  SCR, 
it  will  not  turn  on  and,  thus,  the  voltage  spike  due  to  LI  will  have  no 
effect  on  the  C4  charging  branch.  This  method  does,  unfortunately,  make 
additional  requirements  for  hardware  and  triggering. 

Readjustment  of  the  thyratron  trigger  control  box  settings  is  another 
approach  that  may  be  taken  to  resolve  the  voltage  transient  problem.  The 
firing  mechanism  should  allow  the  voltage  on  the  storage  capacitor  (C4)  to 
exceed  the  inductor  voltage  spike  at  the  time  of  discharge.  This  means 
that  capacitor  Cl  must  be  discharged  at  a  later  time  during  the  half 
cycle.  Thus,  the  new  settings  should  prevent  the  turn-on  of  the  diode 
in  the  discharge  branch  since  it  will  always  be  reverse  biased.  The 
settings  presently  initiate  the  discharge  when  the  source  voltage  reaches 
a  particular  amplitude  during  the  negative  half  cycle. 
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Many  questions  were  left  unanswered  concerning  the  voltage  and  cur¬ 
rent  waveform  anomalies.  They  were  investigated  primarily  to  gain  a 
qualitative  understanding  of  their  underlying  causes.  However,  any 
further  examination  was  considered  to  be  beyond  the  scope  of  this  study. 

Resoi  ince  Determination 

Including  all  of  the  system  charging  inductance  in  the  alternator 
could  possibly  reduce  the  power  source  weight  by  as  much  as  1/3  (Ref  6). 
This  would  remove  the  necessity  for  the  high  leakage  transformer  which 
was  used  in  the  MERADCOM  charging  circuit.  In  an  operational  system  a 
pulse  transformer  between  the  PFN  and  load  cou’d  provide  any  additional 
voltage  multiplication  needed  for  the  discharge  pulse.  Before  the  trans¬ 
former  leakage  inductance  can  be  incorporated  into  the  alternator  its 
internal  inductance  for  the  AC  resonant  mode  should  be  determined.  It 
is  the  combination  of  these  two  inductances  that  must  then  be  designed 
into  the  nigh  reactance  machine.  However,  determining  the  internal 
inductance,  L  ,  turns  out  to  be  a  difficult  task. 

There  is  a  tendency  to  use  classical  steady  state  theory  for  deter¬ 
mining  Lg.  The  initial  step  is  to  find  the  resonant  frequency,  coo> 

The  inductive  and  capacitive  reactances  should  be  equal  at  this  fre¬ 
quency.  Thus, 

*L  =  XC  (7) 

or 

Lc  =  (8) 


where 
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Lc  =  total  charging  inductance 
=  network  charging  capacitance 

Since  L  consists  of  the  internal  alternator  inductance  plus  the  trans- 
c 

former  leakage  inductance,  L  ,  then  L  may  be  written  as 

X  3 


L 

a 


=  L 

c 


-  L 


x 


(9) 


The  difficulty  with  the  above  approach  is  that  classical  steady 
state  resonant  theory  is  not  applicable.  A  popular  circuit  analysis  text 
states,  "Resonance,  by  definition,  is  fundamentally  associated  with  the 
forced  response  since  it  is  defined  in  terms  of  a  purely  resistive  input 
impedance,  a  sinusoidal  steady  state  concept"  (Ref  3:450).  The  response 
of  the  MERADCOM  system  for  half  cycle  charging  includes  both  a  forced 
and  a  transient  component.  This  is  evident  from  the  exponential  terms 
in  equations  (2)  and  (3).  Thus,  the  system  never  reaches  steady  state  in 
the  classical  sense.  This  means  that  the  system  input  impedance  never 
becomes  purely  resistive.  At  the  beginning  of  each  half  cycle  the  source 
"sees"  the  same  input,  a  larr;e  inductor  in  series  with  a  capacitor  charged 
to  a  slightly  negative  value.  If  the  blocking  diode  of  Figure  7  was 
removed,  a  number  of  half  cycles  of  charging  would  be  required  to  get 
the  system  past  the  transient  state.  Only  then  would  the  impedance  be 
resistive  and  the  output  waveforms  sinusoidal.  With  the  present  system 
configuration  only  the  input  is  sinusoidal  and  even  it  is  slightly 
deformed  by  the  turn-on  of  the  diodes.  Therefore,  since  the  source  never 
"sees"  past  the  first  half  cycle,  it  continuously  faces  a  circuit  which 
contains  some  reactance  and  has  a  transient  response.  Thus,  equation  (7) 
is  never  realized. 
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In  most  resonant  circuits  a  number  of  parameters  can  be  varied 
until  resonance  is  detected.  For  example,  changing  either  L,  C,  or  the 
frequency  f  should  yield  the  circuit  response  needed  to  obtain  resonance 
However,  the  frequency  is  usually  chosen  as  the  variable  parameter.  Re¬ 
gardless  of  the  selected  variable,  the  source  terminal  voltage  is  main¬ 
tained  constant.  This  provides  the  standard  from  which  measurements  can 
be  evaluated.  For  example,  resonance  in  the  typical  series  RLC  circuit, 
as  shown  in  Figure  7  with  the  diode  removed,  is  obtained  by  holding  the 
source  terminal  voltage  (point  B)  constant  and  varying  f  until  a  peak 
current  is  reached.  This  occurs  when  the  circuit  becomes  purely  resis¬ 
tive,  that  is,  when  the  inductive  and  capacitive  reactances  cancel  as 
expressed  in  equation  (7).  When  this  happens  a  power  factor  of  unity 
may  be  measured  at  the  source  terminals  and  the  current  will  be  in  phase 
with  the  terminal  voltage. 

The  nature  of  this  research  effort  removes  the  constant  terminal 
voltage  standard  since  part  of  the  charging  inductance  is  inside  the 
alternator.  It  is  iow  desired  to  hold  the  voltage  at  point  A  of  Figure 
7  constant  and  vary  f  until  the  peak  current  is  reached.  Also,  this 
current  should  be  in  phase  with  the  voltage  at  point  A.  La  is  assumed 
to  be  approximately  constant  and  C  is  fixed.  The  difficulty  with  this 
procedure  is  that  Ea,  the  internally  generated  emf  under  load,  is  inac¬ 
cessible.  If  it  could  be  measured  it  would  be  a  simple  task  to  keep  it 
constant  by  varying  the  excitation. 

The  exciter  circuit  for  the  charging  system  is  shown  in  Figure  8. 
The  circuit  enclosed  in  the  dashed  block  rotates  at  the  same  speed  as 
the  rotor.  It  produces  a  flux  which  is  then  conditioned  by  the  armature 
reaction.  The  resulting  flux  then  induces  the  alternator  voltage  Ea 
which  is  approximated  in  the  following  equation. 
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(10) 


Ea  =  4kN<J>f 

where 

k  =  form  factor  to  account  for  non-sinusoidal  aspects 
N  =  number  of  turns  linked  by  the  flux 
<J>  =  peak  amplitude  of  the  resultant  flux 
f  =  impressed  frequency  in  Hz 
Although  equation  (10)  shows  that  is  a  function  of  $  and  f  it  does 
not  infer  that  <f>  is  itself  a  function  of  f.  This  is  the  case,  however, 
since  an  increase  in  f  leads  to  an  increase  in  the  alternator  field 
current  for  constant  exciter  current.  The  field  flux  is  then  boosted 
due  to  this  larger  current.  Thus,  the  total  change  in  results  not 
only  from  the  increase  in  frequency  but  also  from  the  variation  in  flux 
due  to  the  frequency  change. 

From  the  previous  discussion  it  appears  that  there  is  no  direct 

method  for  determining  resonance  and,  hence,  La.  Thus,  an  indirect 

approach  was  taken  which  may  be  used  to  approximate  La> 

At  unity  power  factor  (behind  the  machine  inductance)  Ea  should  not 

be  much  greater  than  the  terminal  voltage  V^,  and  can  be  calculated  if  a 

value  for  L  is  assumed.  Since  E  cannot  be  measured,  this  one  value  is 
a  a 

of  little  help.  However,  if  is  calculated  for  a  number  of  frequencies 
and  is  then  divided  by  the  measured  current  I  for  these  same  frequencies, 
an  impedance  trend  (Z  versus  f)  is  developed.  This  procedure  was  followed 
for  a  number  of  values  for  L  .  The  goal  was  not  to  determine  the  fre- 

o 

quency  at  which  the  impedance  became  purely  resistive  but  rather  the  fre¬ 
quency  at  which  it  became  a  minimum.  As  indicated  earlier,  the  circuit 
is  in  a  transient  mode  and  its  input  impedance  will  always  be  somewhat 
reactive.  The  equation  used  to  calculate  Ea  was 
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(11) 


EM  =  VT/0° 
a  —  1  — 


+  V±  (Ra  +  JVa> 


where 

Ea  =  internally  generated  emf  (volts  rms) 

VT  =  alternator  terminal  voltage  (volts  rms) 

I  =  phase  current  measured  at  the  alternator  terminal  (amps  rms) 

<j>  =  angle  between  and  (degrees) 

0  =  angle  between  and  I  (degrees) 

Ra  =  0.27  =  internal  phase  resistance  of  the  alternator  (ohms) 

=  alternator  frequency  (radians/second) 

La  =  selected  alternator  internal  inductance  (henries) 

The  various  parameters  and  values  which  were  associated  with  equa¬ 
tion  (11)  are  listed  in  Table  I.  All  measurements  were  taken  with  the 
generator  field  current  held  to  a  constant  value  of  8.0  amps.  A  number 
of  values  for  La  were  selected  and  evaluated.  However,  only  three  are 
listed  in  the  table.  The  first  value  selected  was  0.218  mh,  the  speci¬ 
fication  sheet  leakage  inductance.  As  the  frequency  was  varied  the 
impedance  Z  reached  a  minimum  at  2275  RPM.  Although  this  may  indeed  be 
the  frequency  at  which  the  minimum  impedance  occurs  it  is  also  the  fre¬ 
quency  at  which  the  system  started  to  become  slightly  unstable.  It  is 
believed  that  the  60  KW  tests  caused  some  problems  in  the  switching  cir¬ 
cuits.  There  was  insufficient  time  to  try  to  correct  this  problem  and, 
thus,  some  instability  was  noticed  at  the  higher  speeds. 

As  noted  above,  the  choice  of  0.218  mh  for  Lg  resulted  in  a  minimum 
impedance  at  2275  RPM  or  2859  rad/sec.  This  frequency  can  now  be  sub¬ 
stituted  into  equation  (12)  to  determine  the  validity  of  the  L  selection. 

a 


L  c„ 

c  N 


(-Vi1/2 

''2L  '  J 
c 


(12) 
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Data  Used  to  Study  the  Alternator  Response  and  Its  Internal  Reactance 


2350  10.0  40.2  129 


where 


Rc  =  0.445  =  total  charging  circuit  resistance  (ohms)  and  the 
remaining  terms  have  been  previously  defined.  This  equation,  when  solved 
for  Lf ,  yields 


L 

c 


) 


2 


-  (Vi1/2 

CO 

a 


(13) 


Upon  substituting  to  =  2859  rad/sec,  C.7  =  142  pf,  and  R  =  0.445  ohms 
a  N  c 

into  equation  (13),  Lc  is  found  to  be  0.855  mh.  Subtracting  the  trans¬ 
former  inductance  of  0.69  mh  from  this  gives  a  value  for  L  of  0.165  mh. 

a 

This  value  may  now  replace  the  initially  selected  value  of  0.218  mh 

for  L  .  With  this  new  value  E  and  Z  are  calculated  for  varying  fre- 
a  a 

quencies.  As  Table  I  shows,  the  minimum  Z  (both  magnitude  and  angle) 
occurs  at  2275  RPM.  (The  calculated  angle  cannot  be  viewed  quantita¬ 
tively.)  This  same  speed  resulted  for  each  selection  of  L  .  Thus,  the 
alternator  supplies  an  internal  inductance  of  approximately  0.165  mh  at 
a  resonant  frequency  of  455  Hz.  However,  the  data  in  the  last  column 
of  Table  I  implies  that  the  "resonant  bandwidth"  may  be  quite  wide. 

Since  the  alternator  waveforms  are  similar  to  sinusoids,  the  bandwidth, 
as  with  steady  state  theory,  should  depend  upon  the  Q  of  the  circuit. 
This  factor  is  relatively  low  (approximately  5)  and,  therefore,  supports 
the  conclusion  of  a  wide  bandwidth. 

The  value  of  X  ,  determined  from  L  ,  may  now  be  compared  with  the 
values  of  reactance  found  in  the  alternator  specifications  sheet.  For 
this  comparison  it  is  transformed  into  its  per  unit  reactance  value. 

This  is  accomplished  by  the  following  equation 

Xg  =  27rfLa/(Vg/lB)  per  unit  (pu)  (14) 
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where 


f  =  400  =  design  frequency  (Hz) 

Vg  =  240  =  base  voltage  (volts) 

Xg  =  26  =  base  current  (amps) 

Substituting  these  values  and  L  =  0.165  mh  into  equation  (14)  gives  a 

a 

value  of  0.0449  pu  for  X  .  The  only  reactance  value  in  the  specification 
sheets  that  is  in  this  magnitude  range  is  the  armature  leakage  reactance. 
Its  value  is  0.0594  per  unit.  Thus,  the  value  for  appears  to  be 
slightly  lower  than  the  leakage  reactance  value. 

The  alternator  should  respond  differently  for  a  transient  condition 
and  a  steady  state  operation  unless  the  transient  is  of  such  a  nature 
that  it  is  viewed  by  the  machine  as  steady  state.  As  explained  earlier, 
half  cycle  charging  produces  a  transient  response.  So  now  the  point  of 
interest  is  how  this  particular  transient  response  impacts  the  operation 
of  the  alternator.  A  second  use  of  the  measured  data  in  Table  I  may  be 
made  in  order  to  investigate  this  point.  The  data  will  be  applied  to  the 
two-reaction  method  for  computing  the  angle  between  the  phase  current 
Ia  and  the  open  circuit  voltage  (Ref  1:114;  7:204-212).  The  steady 
state  reactance  components  will  be  used  in  conjunction  with  this  data  in 
order  to  develop  the  vector  diagram  drawn  in  Figure  18.  The  two-reaction 
method  assumes  steady  state  operation.  Thus,  if  the  approach  is  still 
applicable  using  the  measured  data  then  it  strongly  suggests  that  the 
alternator  views  the  system  operation  as  steady  state  (nearly  sinusoidal). 
The  following  equations  will  also  be  useful  in  this  development. 
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Figure  18.  Alternator  Vector  Diagram  Using  Steady  State  Reactance  Com¬ 
ponents  (2275  RPM) 


40 


(15) 


I  ,  =  I  sii# 
ad  a 


I  =  I  costJ> 
aq  a 


Xad  *  Xd  -  X£ 


X  =  X  -  x» 
aq  q  t 


(16) 

(17) 

(18) 


where 

I  =  phase  current  (amps) 

^  =  angle  between  and  I  (degrees) 

I  d  =  direct  component  of  the  phase  current  (amps) 

I  =  quadrature  component  of  the  phase  current  (amps) 

X£  =  leakage  reactance  (ohms) 

Xd  =  direct  axis  synchronous  reactance  (ohms) 

Xad  «*  direct  component  of  reactance  (ohms) 

X^  =  quadrature  axis  synchronous  reactance  (ohms) 

Xaq  “  {luadrature  (°r  cross)  component  of  reactance  (ohms) 

The  terminal  voltage  V^,  is  first  drawn  along  the  horizontal.  is 
then  drawn  at  an  angle  9  and  the  1  R  drop  is  laid  off  parallel  to  I  . 

3  3  3 

R  is  the  internal  resistance  of  the  alternator  and  was  measured  to  be 
a 

approximately  0.27  ohms.  From  point  A  the  line  AC  is  drawn  perpendicular 
to  Ia  with  the  magnitude  ^X^.  AB  is  laid  off  along  this  same  line  by 

the  amount  IaX^.  The  angle  between  and  Ia,  is  now  found  by  extend¬ 
ing  a  line  from  the  origin  through  B.  The  magnitude  of  E^  is  determined 

by  segment  CD  which  is  drawn  perpendicular  to  this  line.  This  value 

should  closely  approximate  the  measured  open  circuit  voltage  if  the 
alternator  may  be  viewed  as  operating  in  a  steady  state  mode. 

The  direct  and  quadrature  components  of  the  phase  current  may  now 
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be  used  to  find  the  internally  generated  voltage,  Ea.  This  is  begun  by 
first  extending  a  line  from  point  D  to  point  F.  It  is  parallel  to  CD  and 
has  a  magnitude  Iaq^aq*  The  vector  IacjXa£j  is  then  drawn  perpendicular  to 
this  line.  It  begins  at  point  F  and  should  terminate  at  some  point  G  on 
line  segment  AB.  The  vector  for  the  internally  generated  voltage  Eg  may 
then  be  drawn  from  the  origin  to  G.  The  final  vector  of  interest  connects 
points  A  and  G  and  should  have  a  value  I  Xp . 

3. 

Many  cystem  transients  are  of  such  a  nature  that  they  induce  tran¬ 
sient  currents  inside  the  alternator.  For  this  case  the  synchronous 
reactance  values  cannot  be  used  in  the  development  of  the  vector  diagram 
in  Figure  18.  A  use  of  these  values  would  cause  the  diagram  to  be  in 
error.  Most  notable,  E^  would  be  incorrect  in  magnitude  and  direction 
and,  hence,  would  not  equal  its  measured  value. 

Under  certain  circumstances  the  operation  of  the  alternator  may  be 
viewed  as  steady  state  even  though  it  feeds  a  circuit  experiencing  a 
transient  response.  Apparently,  this  is  what  is  happening  during  AC 
resonant  charging  and,  thus,  the  use  of  the  synchronous  reactances  in 
the  development  of  the  alternator  vector  diagram  is  permissible. 

There  are  a  number  of  reasons  the  alternator  may  be  considered  to 
be  operating  in  the  steady  state.  At  the  beginning  of  each  charging  half 
cycle  the  circuit  connected  to  the  alternator  terminals  is  changed.  Just 
before  this  half  cycle  begins  the  source  is  completing  the  build-up  of 
peak  charge  on  one  capacitor.  At  the  completion  of  this  build-up  the 
charged  capacitor  circuit  is  replaced  with  a  circuit  containing  a  capa¬ 
citor  holding  a  slightly  negative  charge.  This  new  circuit  results  in  a 
mathematically  transient  response  (see  equation  (2)).  However,  the 
response  is  so  gradual  (primarily  due  to  the  large  charging  inductance) 
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that  the  alternator  can  barely  distinguish  it  from  400+  Hz  operation.  As 
shown  in  Figure  12,  it  contains  no  current  discontinuities  at  the  half 
cycle  end  points  or  during  the  charging  of  the  capacitor.  The  smooth 
transition  at  the  end  points  is  primarily  due  to  the  initial  condition 
requirements.  These  are  specified  in  Appendix  A  which  shows  that  the 
current  must  be  zero  at  both  the  beginning  and  end  of  a  half  cycle.  All 
of  these  characteristics  combine  to  produce  a  repetitive  waveform  at  the 
alternator  terminals  and  an  armature  reaction  that  is  fixed  both  in  magni¬ 
tude  and  direction  with  respect  to  the  poles.  The  fixed  nature  of  the 
armature  reaction  prevents  the  build-up  of  the  voltages  which  cause 
transient  currents  inside  the  machine  (Ref  5:286-290).  Thus,  synchronous 
reactances  may  be  used  in  developing  a  two-reaction  alternator  vector 
diagram  as  shown  in  Figure  18.  The  vectors  drawn  in  this  diagram  were 
based  upon  Table  I  measurements  taken  at  2275  RPM  and  the  synchronous 
alternator  reactances  listed  in  Appendix  D.  The  diagram  results  in  an 
open  circuit  voltage  of  215  volts.  The  terminal  point  of  the  Ia{jXad 
vector  is  slightly  past  line  AB.  The  distance  from  this  point  to  A  is 
measured  to  be  approximately  24.6  volts  which  results  in  a  value  for 
of  0.064  per  unit.  This  per  unit  value  is  just  slightly  larger  than  the 
0.0594  leakage  reactance  value  specified  for  the  alternator.  Similarly, 
the  graphical  value  of  Ej  is  6  volts  less  than  the  measured  value  of  221 
volts.  These  differences  are  well  within  the  realm  of  error  in  measure¬ 
ment.  The  results  indicate  that  although  half  cycle  charging  is  a 
transitory  condition  it  allows  the  alternator  to  experience  a  mode  of 
operation  that  very  closely  resembles  steady  state. 
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IV.  Conclusions  and  Recommendations 


There  were  a  number  of  objectives  associated  with  this  study.  One 
of  the  most  important  was  to  determine  the  impact  on  an  alternator  when 
it  operates  in  the  AC  resonant  charging  mode.  Its  output  waveforms  were 
to  be  analyzed  for  degradation  while  /arying  the  load  and  frequency.  The 
final  objective  was  to  determine,  if  possible,  when  the  system  was  at 
resonance.  This  was  equivalent  to  determining  the  internal  reactance  of 
the  machine. 

Conclusions 

The  system  was  successfully  tested  for  15,  30,  and  45  KW.  However, 
only  a  few  measurements  were  taken  at  the  60  KW  level  due  to  the  develop¬ 
ment  of  trigger  circuit  problems.  These  load  variations  did  not  cause 
any  noticeable  degradation  of  the  alternator  operation.  A  similar  con¬ 
clusion  was  reached  when  the  frequency  was  varied  from  approximately  380 
to  475  hertz.  This,  to  some  extent,  results  from  a  low  circuit  Q  of 
approximately  5. 

Although  a  number  of  anomalies  were  detected  in  the  charging  wave¬ 
forms,  the  waveshapes  very  closely  conformed  to  the  theoretical  predic¬ 
tions.  The  anomalies  were  found  to  be  peculiar  to  the  hardware  and  con¬ 
figuration  of  the  system  under  investigation. 

During  the  latter  tests,  system  noise  tended  to  aggravate  the 
trigger  circuits  and  caused  pre-firing  of  certain  thyratrons  that  must 
have  been  adversely  affected  during  the  60  KW  tests.  Reducing  the  gas 
pressure  inside  the  thyratron  tubes  could  have  possibly  alleviated  this 
condition. 


i 
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No  straightforward  method  was  developed  for  determining  precise 
resonance.  The  term  resonance  as  used  in  this  study  context  is  some¬ 
what  of  a  misnomer  since  classical  steady  state  resonant  theory  is  not 
really  applicable.  Classical  theory  addresses  only  the  forced  response. 
This  study  showed  that  the  transient  response  cannot  be  excluded. 

Classical  theory  also  seeks  a  maximum  amplitude  response  when  all 
reactances  cancel  and  the  system  is  purely  resistive.  This  study  con¬ 
cludes  that  the  half  cycle  charging  circuit  has  a  transient  response. 

The  reactances  do  not  completely  cancel  and,  thus,  the  input  impedance 
is  never  purely  resistive.  Even  so,  the  alternator  armature  reaction 
is  approximately  fixed,  both  in  magnitude  and  direction,  with  respect 
to  the  poles.  This  implies  that  the  alternator  is  in  a  near  steady  state 
mode  of  operation.  This  latter  point  was  confirmed  by  the  use  of  measured 
data  and  synchronous  reactances  in  a  two-reaction  vector  analysis  for  the 
alternator.  This  same  data  was  used  to  determine  the  internal  reactance 
of  the  alternator  during  resonant  charging.  It  was  found  to  be  slightly 
less  than  the  armature  leakage  reactance. 

Recommendations 

The  literature  search  for  this  study  led  to  the  conclusion  that  there 
is  only  scant  material  on  charging  systems.  Thus,  an  investigation  of 
methods  other  than  AC  or  DC  resonant  charging  would  contribute  to  filling 
this  void.  One  method,  in  particular,  which  has  caught  the  interest  of 
the  Army  (MERADCOM)  is  network  charging.  This  method  limits  the  peak 
current  requirements  and  is  not  expected  to  cause  as  negative  an  impact 
on  the  alternator  as  does  DC  resonant  charging.  Further,  its  controls 
and  triggering  should  be  simpler  than  those  needed  for  AC  resonant 
charging.  It  needs  to  be  investigated  and  compared  with  the  methods 
already  studied. 
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The  AC  pulser  is  a  very  basic  prototype.  The  MERADCOM  engineers 
were  aware  of  many  of  its  shortcomings  before  its  delivery.  If  these 
design  problems  could  be  corrected  this  would  permit  more  reliable  and 
expanded  testing.  In  addition  to  redesigning  the  trigger  controls,  con¬ 
sideration  would  have  to  be  given  to  wiring,  grounding,  noise  suppres¬ 
sion,  and  the  general  physical  layout  of  the  circuit. 
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Appendix  A 


Theoretical  Development  of  the 
AC  Resonant  Charging  Equations 
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.  .  is  taken  from  Ref  2.  The  differential 

The  following  analysis  is  tak 

f  Fiaure  7,  in  terms  of  the  network  charge 
equation  for  the  circuit  of  Figure  , 


V  is 


T  d  qN  .  R  +  -A-q  =  E  cos  (tot  +  4>) 

Vtt  +  Rc  it  +  a 


The  s[«aay  «*.  oS  <-»  -  «“  ,m 


qi(t)  -  Q„  EJ“»e 


where  QN  Is  found  by 


substituting  <A-2>  into  (A-l).  Thus, 


CN  Ea£ 


,N  '  1  -  Wa2  *  3  Wa 


coulombs 


transient  integral  for  (A-l)  is 


q^(t)  =  A  £pt  coulombs 


where 


A  is  a 


complex  constant  of  integration 


K  /  1  ,  C  y 

P  =  _a  +  jia=-2iT+iV  ITc^  2Lc 

complete  solution  for  (A-l)  is  then 


j<D  t  ,  7  _Pt 

,„(t)  •  <.!«)  +  12<‘>  *  %  E  ‘  +  *  " 


Thus,  equation  (A-3)  may  be  written  as 


(A.-12) 


qN(t)  =  Q^cosCa^t  +  <$)  -  9)  +  A  e  cos(o)at  +  i})) 

Expanding  cos(<bat  +  <p)  and  recognizing  that  =  AcosiJ;  and  A2  =  Asin^ 
(A-12)  becomes 

qN(t)  =  QNsin(o)at  +  <p)  +  £  [A^cosa^t  -  A2sincoat] 


A^  and  A2  in  this  equation  can  now  be  replaced  with  the  right  side  of 
equations  (A-7)  and  (A-8) .  Making  these  changes  and  substituting  Q^sin<(> 
for  and  Q^cos^  for  Q2  yields  the  equation  for  the  network  charge 


qN(t)  =  Qn(1  -  £  a£)sin(uat  +  <j>)  +  e-at{qN(0)cosa)at  + 


i  (0)  aqN(0) 

(—5“  +  ~ - ^sittNsinw  t)  (A-13) 

a  a  a 


The  condition  for  repeating  transients  is,  for  half  cycle  charging, 
i c (it)  +  ic(0)  =  0.  The  current  is  now  obtained  by  differentiating  equa¬ 
tion  (A-13)  and  is 
E 

ic(t)  ■  ip  (1  “  e_at)cos(oJat  +  4>) 
c 

E 

+  e  at{i  (0)cosu  t  -  cu  q,T(0)sino>  t  +  — sinw  t  cos  <j>) 

c  a  a^N'  a  2L  oi  a 

c  a 

E 

+  e_a  {2lV  sin*~  ic(0)  "  a<5N(0)}  sina)at  (A_U) 

a  ^c  a 

E 

where  has  been  substituted  for  Q^. 

c  a 

At  <i)  t  =  it  equation  (A-1A)  becomes 

3. 


f 

1 


I 
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i  (tt)  = 
c 


E  -  15.  _  15. 

-2.  (1  -  e  «a)cos<t>  "  E  “a 
c 


-  V°> 


Solving  for  i  (°)  yields 


1(0)  =  cos* 
c  Rc 


(A- 15) 


The  quality  factor  Q  is  defined  as 


Lc“a 

O'  — 


R  c„“ 

c  c  N  a 


(A- 16) 


substituting  1(0)  and  Q  into  (*-13)  a  revised  equation  for  the  net- 


By  ouuow** - w  c 

work  charge  is  obtained 


6)  t 
a 


qN(t) 


q  cNEa  (1  -  e*  2Q)sin(wat  +  *) 


CO  t 

a 


C„E 


4.  J^,._(o))sinco  t} 


+  e  2QlqN(0)coscoat  +  (QC^cos*  -  ~r^ln*  +  ^ 


a 

(A-17) 


ttu  e,*,^  IS  w  P—  *.  -»««. 


OJ  t 
a 


i  <t)  .  QCaVa  «  -  E'  2Q>oos(“at  +‘ 


(0  t 
a 


+  e~  2Q  {QCNEa<oacos4>cosa)at 


,qN 


(0) 


"  lC//alCE  ”  4Q‘ 
n  a 


iin*) 


qM(o)<o  i 

_  J? — — 2-lsincoat} 

4Q 


(A-18) 


This 


expression  may  then  be  restructured  to  give 
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t  _ 

l 

E  t 

i  (t)  =  -rf—  cos(co  t  +  4>)  +  i  (O)cosco  t 
c  L  L  3  c  a 

c 

I 

* 

E 

+  [2L  co  COS0  "  “a  qN(0>Jsin“at 

C  <3 

(A-22) 

Thus,  for  half  cycle  charging 

I 

*cM  "  -  coscj)  +  i<0)] 

a  c 

(A-23) 

( 

As  stated  earlier,  the  condition  for  repeating  transients  is 

ic(TT)  - 

-  ic(0).  Therefore,  the  only  solution  for  (A-23)  is 

f 

cos$  *  0  or  =  ±  -j 

and  so,  from  (A-15) 

i  (0)  =  0 
c 

1 

Equation  (A-21)  may  be  evaluated  at  the  discharge  point 

to  determine  the  resulting  stepup  voltage  ratio 

V  = 

r 

vn(tO  _  jl  qN(o>  _  — 

.  -  [Q(1  _  e  2Q)sin<f>  +  c  _  e  2Q] 
a  N  a 

(A-24) 

\ 

If  tt/2Q  <<  1  then  the  exponentials  in  the  above  equation  may 

be  approxi- 

•*' 

mated  using 

e"  ^  *  1  -  ^ 

and 
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1 

1  . 

% 

•• 

/  * 

■  :  i 

1  -  e 


IT 

2Q 


—  (1  _  JL) 
2Q  4q' 


Substitution  of  these  expressions  into  (A-24)  yields 


vn(tt} 


-  tf  (1 


4Q 


•)sin<f  + 


qN(0) 

CNEa 


(1  -^)] 


(A-25) 


Then  for  the  case  where  there  is  no  initial  charge  on  the  network  and 
the  losses  are  negligible,  (A-25)  becomes 


-f—  *  -  ?  sin* 


(A-26) 


and,  thus,  for  maximum  network  voltage 


<t>  =  - 


TT 

2 


(A-27) 
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PROGRAM  EQUATE 


100-C  THIS  PROGRAM  INCLUDES  THE  EQUATIONS  FOR  THE  CHARGING  VOLTAGE 
U0=C  AND  CURRENT  AT  RESONANCE.  THE  SOURCE  FOR  THE  EQUATIONS  IS  THE 
120=C  PULSE  GENERATOR  TEXT  BY  GLASOE  AND  LE3ACQZ.  FOR  COMPARISON 
130=C  PURPOSES  THE  SOURCE  AMD  INDUCTOR  VOLTAGE  EQUATIONS  HAVE  ALSO 
140=C  BEEN  INCLUDED.  THE  TIME  STEP  FOR  THE  EQUATIONS  HAS  BEEN  SET  TO 
150=C  IO.OE-6.  THE  EQUATIONS  ARE  SOLVED  AT  THESE  TIME  POINTS  FOR 
160=C  A  HALF  PERIOD.  RATHER  THAN  DISPLAYING  THE  OUTPUT  ON  THE 
I70=C  SCREEN,  IT  IS  STORED  ON  TAPE8. 

130=C 

190=C 

200=C 

210=  PROGRAM  GLASOE (INPUT, OUTPUT, TAPE5=INPUT,TAPC6=OUTPUT, TAPES) 

220=  REAL  LC , IC , ICMAX 

230=  REWIND  8 

240=  PI=3. 141592654 

250=  RC=0.445 

260=  LC=0.353E-3 

270=  CN=147.17F-6 

280=  EA=150.0 

290=  VCZERO=-5.4 

300=  Plll=-P  1/2.0 

310=  A=RC/(2.0*LC) 

320=  OMEGAO= 1 . 0/ SQRT ( LC*CN ) 

330=  OMEGAA=SQRT(OMEGAO**2-A**2) 

340=  RPM=OMECAA/(0.4*PI) 

350=  Q=OMEGAA*LC/RC 

360=  QCZERO=CN*VCZERO 

370=  F I NT IM=  PI/ OMEGAA 

380=  T-0.0 

390=  WRITE(6 ,3) 

400=  3  FORMAT( 6  X , 4HTIME , 9X , 2HI C , 1 2X , 2HVC , 1 OX , 2HVL, 10X , 6H0I1EGAA) 

410=  1  EX=EXP(-0MEGAA*T/(2.0*Q)) 

420=  IC=Q*CN*EA*0MEGAA*((1.0-EX)*C0S(0'1EGAA*T+PHI) 

430=  >/+EX*(COS(?HI)*COS(OMEGAA*T) 

440=  /*-(  QCZERO/  (  Q*CN*EA)*(  1 .0- 1 .0/  (  4 .0*Q**2)  ) 

450=  #~1 .0/ (4 ,0*Q**2)*SIN(PHI ) )*SIN(OMEGAA*T) ) ) 

460=  VC=EA*0*( ( 1 ,0-EX)*SIM(OMECAA*T+PHI)-t-EX* 

470=  H QCZERO/ (Q*CH*EA)*COS(OMEGAA*T)+ 

480=  #(C08(PHI)-1. 0/(2. 0*Q)*SIN(PHI)+1. 0/(2. 0*Q**2)* 

490=  #QCZF,RO/(CN*EA))*SIN(OMEGAA*T))) 

500=  V=EA*C0S(0MEGAA*T+P!1I ) 

510=  VL=V-RC*IC-VC 
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520=* 
530=0 
540=  2 
550»  4 
560=  5 
570= 
530= 
590= 
600= 


WRITE (G ,2)T, IC,VC ,VL,OMKGAA 
WRITE (6 , 2)T, IC,VC ,VL,OMEGAA 

F0RMAT(1X,E12.5,1X,E12.5,1X,E12.5,1X,E12.5,1X,EI2.5) 

ICHAX=IC 

T=T+10.0E-6 

IF(T.LE.FINTUt)GO  TO  1 

REWIND  8 

STOP 

END 


r*- 
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Appendix  C 

Operating  Instructions  and  Listing  for  Program  RKPOUR 
Program  RLCCHG 
Program  RLCDIS 
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Operating  Instructions  for  RKFOUR 

Programmer:  B.  D.  Weathers 

Purpose:  RKFOUR  is  an  integration  routine  designed  to  provide  the 

equivalent  of  an  analog  computer  integrator  for  the  solution 
of  normal-form  differential  equations. 

Features 

-  An  arbitrary  number  of  integrators  can  be  used.  The  number  used  is 

determined  by  the  size  of  the  storage  arrays  set  in  the  calling 
program. 

-  All  variables  that  affect  the  operation  of  RKFOUR  appear  in  the 

calling  statement;  multiple  uses  of  the  program  are  possible,  but 
care  must  be  used  to  prevent  time  desynchronization. 

-  The  independent  variable  step  size  is  determined  automatically,  but 

if  desired  this  program  can  be  made  to  run  with  a  fixed  step  size. 

-  The  adaptive  step  size  feature,  if  used,  will  insure  little  trouble 

with  numerical  instability. 

-  The  limits  on  the  adaptive  step  size  may  be  set  in  the  calling 

program,  or  default  values  may  be  used. 

-  The  error  limits  that  determine  when  the  step  size  should  be  changed 

are  program  parameters,  and  hence  may  be  externally  controlled. 

-  An  adjustable  delay  is  available  to  prevent  eratic  step  size 

changing.  Once  the  criteria  is  satisfied  step  size  increases  by 
the  appropriate  factor  at  each  time. 

-  Step  sizes  may  be  changed  externally,  but  may  not  be  changed  in  the 

middle  of  an  iteration,  that  is  when  KEEP  is  equal  to  0. 
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A  solution  point  counter  is  contained  in  the  routine  to  enable 
termination  if  appropriate  independent  variable  criteria  is  not 
available. 

Suggestions  for  Use 

Before  using  this  subroutine  for  the  first  time,  a  careful  examina¬ 
tion  of  the  flow  chart  should  be  made.  Also  a  thorough  examination  of 
the  program  parameter  should  be  carried  out  to  insure  that  the  default 
values  of  the  parameters  will  be  satisfactory,  if  used. 

You  must  set  all  non- zero  initial  integrator  outputs.  All  non-unity 
integrator  gains  must  also  be  set  in  section  2. 

The  value  for  H(l)  -  Step  Size  -  should  be  set  as  close  to  the  pre¬ 
dicted  sampling  interval  as  possible.  This  will  insure  a  more  efficient 
operation  for  the  first  pass. 

If  running  in  "Open  Shop"  and  the  20  second  time  limit  is  just  a 
"bit  too  short,"  add  "RUN=N0  CHECK"  to  the  control  card.  Do  this  only 
after  your  name,  after  you  are  SURE  the  program  is  free  of  programming 
errors.  On  a  large  job  this  can  add  15%  to  the  available  computation 
time,  by  eliminating  several  compiler  tests. 

CALL  RKFOUR  (NINT,  Y,  X,  GA,  PE,  XX,  H,  ZZ1,  ZZ2,  ZZ3,  TIME,  POTS,  KEEP, 
IFL) 

NINT  -  Number  of  integrators.  Storage  for  the  integrators  must  be 
defined  in  the  calling  program.  There  is  no  maximum,  but 
the  example  program  is  set  for  up  to  20. 

Y  -  An  array  that  contains  the  integrator  output.  Default 
initial  values  are  zero. 

X  -  An  array  that  contains  the  integrator  inputs. 
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GA  -  An  array  that  contains  constant  gains  associated  with  each 


integrator.  Default  values  are  1.0. 

PE  -  A  working  array  used  to  store  absolute  peak  values  of  the 

integrator  outputs.  This  array  is  set  to  zero  during  initial 
entry  to  RKFOUR. 

XX  -  An  array  used  in  two  ways.  If  I  denotes  the  Ith  integrator, 
then  XX(1,1)  stores  the  past  value  of  the  integrator  output 
and  XX(2,I)  stores  the  past  value  of  the  integrator  input. 
During  the  integration  cycle  the  XX(3,I)  and  XX(4,I)  positions 
are  used  to  store  the  derivative  values  needed  by  the  integra¬ 
tion  algorithm.  At  the  end  of  a  cycle,  XX(3,I)  holds  the  new 
integrator  output  and  XX(4,I)  holds  the  new  integrator  input. 
It  is  possible,  with  the  values,  to  implement  a  third  order 
interpretation  for  the  values  of  the  integrator  outputs 
between  steps.  This  would  serve  to  offset  the  problem  of 
obtaining  values  at  discrete  time  points  and  would  be  more 
efficient  than  step  size  manipulation. 

H  Array  -  The  integration  control  parameters. 

H  (1)  -  The  current  independent  variable  step  size. 

H  (2)  -  Minimum  permissible  step  size  (Default  value  is  0.0). 

H  (3)  -  Maximum  step  size  (Default  value  is  1.0E  +  8). 

H  (4)  -  Factor  used  to  modify  the  step  size— must  be  greater 
than  1.0  (Default  value  is  SQRT  of  2=1.414214). 

H  (5)  -  Proportional  to  the  lower  error  bound  (Default  value 
is  1.0E  -  4). 

H  (6)  -  Proportional  upper  error  bound  (Default  value  is 
1.0E  -  02). 
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H  (7)  -  Adaptive  step  size  control  (Set  to  -1  for  constant 
step  size  operation)  (Default  value  of  0.0  causes 
adaptive  step  size  operation). 

H  (8)  -  Delay  in  increasing  the  step  size  after  error 
criteria  is  met.  (Default  value  is  3.0). 

H  (9)  -  Last  value  of  the  independent  variable. 

ZZ1,  ZZ2,  ZZ3  -  Working  arrays  used  to  store  the  constants 

needed  for  the  integration  routine.  These  arrays  are 

TIME  -  The  independent  variable  set  within  RKFOUR.  The 
default  initial  value  is  zero. 

PNTS  -  The  counter  of  permanent  solution  extensions.  Use 

this  variable  for  termination  when  other  termination 
options  fail. 

KEEP  -  This  variable  is  set  to  ONE  only  when  the  data  in  the 
derivative  evaluation  section  is  permanent.  It  is 
zero  at  all  other  times.  This  variable  should  be  used 
as  a  control  upon  memory  dependent  functions . 

IFL  -  The  IFLag  array. 

IFL(l)  -  Set  to  -1  at  the  first  entry  to  RKFOUR,  and  is  0  at 

the  first  exit.  On  subsequent  exits  it  has  a  value  +1. 

IFL(2)  -  Has  the  value  -1  to  +3  depending  upon  where  the  inte¬ 
gration  routine  is  in  its  operation  cycle. 

IFL(3)  -  This  flag  is  used  to  control  the  step  size  from  the 
derivative  evaluation  section  of  the  calling  program. 

It  is  set  to  -1  for  step  size  reduction,  and  +1  for 
an  increase  in  the  step  size. 

IFL(4)  -  This  flag  is  -1  when  initializing  calculations  are  to 
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be  made.  It  is  0  when  only  the  derivative  evalua¬ 
tion  is  required,  and  is  +1  when  the  solution  has 
been  successfully  extended  over  one  step  of  the 
independent  variable. 

IFL(5)  -  This  flag  is  set  to  +1  immediately  after  IFL(4)  is 
-1  for  the  first  time.  It  is  intended  as  a  signal 
for  other  routines  that  may  be  used  that  initializa¬ 
tion  is  underway. 

IFL(6)  -  This  flag  is  used  for  the  temporary  storage  of  the 


variable  KEEP. 


The  Basic  RKFOUR  Program 


DIMENSION  Y(20)4,  X(20)4>  GA(20)4,  PE(20),  XX(4,20),  H(9) 
DIMENSION  ZZ1(20) ,  ZZ2(20),  ZZ3(20),  IFL(7) 

If  Other  Dimensioned  Variables  are  required  the  necessary  DIMENSION 
statements  go  here. 


SECTION  1  INSERTIONS  GO  HERE 

NRUN  =  (Number  of  runs  to  be  made) 

DO  6  IRUN=1,NRUN 
ISTOP=0 
KEEP=1 
IFL(1)=-1 

NINT= (Number  of  integrators  required) 

1  IFL(2)=-1 

2  CALL  RKF0UR(NINT,Y,X,GA,PE,XX,H,ZZ1,ZZ2,ZZ3,TIME,PNTS,KEEP,IFL) 
IF  (IFL[4])3,4,5 

3  INT=1 
TIME1=-1. 256789 
NC=1 

NCC=1  (Portion  of  Points  released  to  Section  5) 

FINTIM=  (Maximum  Value  for  TIME  in  this  run) 

PTMAX=  (Maximum  Number  of  solution  points  for  this  run) 


SECTION  2  INSERTIONS  GO  HERE 

4  IF(TIME  -  TIME1)7,8,9 
7  TIME1=TIME 


SECTION  3  INSERTIONS  GO  HERE 
8  CONTINUE 


SECTION  4  INSERTIONS  GO  HERE  (Exclude  everything  nonessential) 

GO  TO  2 
5  IPR=0 
NC=NC-1 
IF  (NC)9, 9, 10 
9  NC-NCC 
IPR=1 

10  IF (TIME . GE . FINTIM . OR . PNTS . GE . PTMAX)  IST0P=1 
IF  ( [IPR+ISTOP ] .EO.O)  GO  TO  11 


SECTION  5  INSERTIONS  GO  HERE 


11  INIT=0 

IF  (ISTOP)l,l,6 
6  CONTINUE 

WRITE  (6,100) 

100  FORMAT  (1H1,  IX) 

STOP 

END 

Basic  RKFOUR  System 

There  are  a  number  of  variables  in  the  RKFOUR  System  whose  value 
depends  upon  the  problem  undertaken. 

NRUN  -  The  number  of  runs  to  be  made. 

NINT  -  The  number  of  Integrators  to  be  used  in  the  run  (limited 
to  20  by  storage  allocation  in  the  stock  system). 

NCC  -  The  portion  of  solution  points  that  are  released  to  the 
output  section  (Section  5)  NCC  =  1  delivers  every  point, 
NCC  =  2,  every  second  point,  etc. 

FINTIM  -  The  maximum  value  that  the  independent  variable  will 

attain.  The  program  will  terminate  upon  the  first  solu¬ 
tion  point  to  meet  or  exceed  this  value. 

PTMAX  -  The  maximum  number  of  solution  points  in  any  one  run. 

IPNT  -  The  point  counter  used  in  plot  routines.  It  may  be  set 
to  0  in  one  of  two  locations.  If  summary  curves  are 
desired,  set  to  0,  before  the  IRUN  DO  statement.  If  no 
curves  are  desired,  set  to  0  after  the  IRUN  DO  statement. 

INIT  -  This  variable  is  1  the  first  interation  step,  and  0 

thereafter.  This  can  be  used  as  a  memory  flag  for  memory 
dependent  functions. 

The  Basic  RKFOUR  System  contains  five  sections  where  the  particular 
problem  statements  may  be  inserted. 


67 


1.  This  section  is  used  to  initialize  those  values  which  are  not 
run  dependent,  or  whose  calculation  in  block  form  is  more 
efficient.  In  this  case  they  may  be  indexed  by  IRUN  for  use 
in  the  rest  of  the  program. 

2.  This  section  is  used  to: 

a.  Set  non-zero  integrators  initial  conditions 

b.  Set  non-unity  integrator  gains 

c.  Override  any  default  values  set  by  RKFOUR  subroutine 

d.  Carry  out  calculations  that  are  run  dependent  and  need  be 
carried  out  only  once  per  run 

3.  Thx.  T’ction  is  used  to  carry  out  all  calculations  of  the  func¬ 
tion^  -'j  the  independent  variable. 

4.  This  section  contains  the  derivative  evaluation  equations.  It 
is  executed  four  times  for  each  point  and  should  be  carefully 
structured.  For  example: 

X(1)  =  (A+B) /(IH-E)  +  U(1)*(G+H)  would  be  efficiently  written  as 
X(l)  *>  AB  +  Y(l)  *  GH 
where  AB  =  (A+B) /(DfE) 

GH  «  G  +  H 

and  are  evaluated  in  section  1,2,  or  3  as  is  appropriate. 

5.  This  section  is  used  to  call  the  output  routines.  If  NCC  is 
not  equal  to  1  it  is  not  entered  on  every  pass,  and  this  should 
be  considered  when  structuring  the  problem. 
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SUBROUTINE  RKFOUR 


100=  SUBROUTINE  RKFOUR(N,Y,X,CAXN,PEAK.,XX,U,Cl ,C2,C3 ,TIME,PNTS,K 

EEP.FL) 

110=  DIMENSION  Y(20),X(2G),GAIN(20),PEAK(20),XX(4,20), 

120=  $H(9) ,Cl(20) ,C2(20) ,03(20) ,A(8) 

130=  INTECER  FL(7) 

140=  DATA  A(1)/0.1/,A(2)/0.0/,A(3)/1.0E8/,A(4)/1. 414214/, A(5)/l. 

03E-4/ 

150=  DATA  A(6)/1.0E-2/,A(7)/0.0/,A(8)/3.0/ 

160=C 

170=  IF(FL(1))  1,3,20 

180=1  DO  2  1=1, H 

190=  Y(1)=0.0 

200=  PEAK(I)=0.0 

210=2  GAIN( I)=l .0 

220=  TIME-0.0 

230=  PNTS-0.0 

240=  FL(4)=-1 

250=  FL(5)=I 

260=  FL( 1 )=0 

270=  DO  9  I»  1,8 

280=9  H(I)=A(1) 

290=  RETURN 

300-3  Q=H( 1 ) 

310=  FL(4)»1 

320=  FL(7)=H(8) 

330=4  H(1)=Q 

340=  FL(3)*0 

350=  DO  5  1=1, N 

360=  C2(I)=GAIN( I)*Q 

370=  C1(I)=0.5*C2(I) 

330=5  C3(I)=C2(I)/6.0 

390=  IF(FL(1))  6,6,10 

400-10  IF(FL(2))  31,33,33 

410=6  FL( 1 )=1 

420=7  PMTS-PHTS+1.0 

430=  DO  8  1=1, N 

440=  XX(3 ,I)=Y(I) 

450=8  XX(4,I)=X(I) 

460=  FL(4)-1 

470=  H(9)=TIME 

480=  RETURN 

490=20  IF(FL(3))  25,30,25 

500=25  IF(FL(2)-3)  26,60,60 

510-26  IF(FL(3) )  21,30,22 

520=21  q»H(l)/H(4) 

530=  IF(Q“H(2))  23,4,4 

540=23  Q-H(2) 

550=  IF(Q-H(1))  4,30,4 

560=22  Q=H(1)*H(4) 

570=  IF(0-H(3)>  4,4,24 

580=24  Q=H(3) 


590=  IF(Q-llO))  4,30,4 

600=30  FL(5)=0 

610=  IF(FL(2) )  31,40,45 

620=31  FL(6)=KEEP 

630=  KEEP“0 

640=  DO  32  I*1,N 

650=  XX(1,I)=XX(3,I) 

660=32  XX(2,I)=XX(4,I) 

670=33  FL(2)=0 

680=  TIME-  H(9)+0.5*H(1> 

690=  DO  34  1=1, N 

700=34  Y(I)=XXa,I)+C!<I>*XX(2,I> 

710=35  KEEP  =  0 

720=36  FL(4)=0 

730=  RETURN 

740=40  FL(2)=1 

750=  DO  41  1=1, N 

760=  XX(3,I)=X(I) 

770=41  Y(I)=XX(1,I)+C1(I)*X(I) 

780=  GO  TO  35 

790=45  IF(FL(2)-2>  46,48,60 

800=46  FL(2)=2 

810=  TIf  !E=11  ( 9)+H  ( 1 ) 

820=  DO  47  1=1, N 

830=  XX(4,I)=X(I) 

840=47  Y(I)=XX(1,I)+C2(I)*X(I) 

850=  GO  TO  35 

860=48  FL(2)=3 

870=  IF(H( 7 ) )  49,52,52 

880=49  KEEP=FL(6) 

890=50  DO  51  1=1, N 

900=51  Y(I)=XX(l,I)+C3a)*(XX(2,I)+X(I)+Z.*(XX(3,I)+XX<4,I))) 

910=  GO  TO  36 

920=52  ERR=0.0 

930=  DO  56  1=1, N 

940=  ESTERR-  C2(1)*ABS(XX(2,I)+X(I)-2.*XX(4,I)) 

950=  IF(PEAKU))  53,53,54 

960=54  ESTERR=ESTERR/PEAK(I) 

970=53  IF (ERR- ESTERR)  55,56,56 

980=55  ERK=ESTERR 

990=56  CONTINUE 


1000=  IF(ERR-H(6))  57,57,66 

1010=57  IF(ERR-H(5))  58,59,59 

1020=58  FL(7)=FL(7)-1 

1030=  IF(FL(7) )  65,65,49 

1040=65  FL(3)=1 

1050=  GO  TO  49 

1060=59  FL(7)=H(8) 

1070=  GO  TO  49 

1080=66  FL(7)=H(8) 

1090=  GO  TO  21 

1100=60  IF(H(7))  7,61,61 

1110=61  DO  63  1=1, N 

1120=  ERR=ABS(Y(I)) 

1130=  IF(ERR~PEAK(I))  63  ,63,62 

1140=62  PEAK(I)*ERR 

1150=63  CONTINUE 

1160=  GO  TO  7 

1170=  END 


PROGRAM  RLCCHG 


100=C  THIS  PROGRAM  USES  THE  RKFOUR  SUBROUTINE  TO  INTEGRATE  THE 
110=C  DIFFERENTIAL  EQUATIONS  OF  AN  RLC  CIRCUIT  WITH  A  SINUSOIDAL 
120=C  INPUT  AND  AN  INITIAL  VOLTAGE  ON  THE  CAPACITOR.  THE  CHARGING 
130-C  CURRENT,  CHARGING  VOLTAGE,  AND  INDUCTOR  VOLTAGE  ARE  EACH 
140=C  CALCULATED  AND  WRITTEN  TO  TAPES  WHERE  THEY  MAY  BE  EDITED 
150=C  AND  VIEWED.  THE  INTEGRATION  TIME  IS  OVER  A  HALF  CYCLE. 

160=C 

170=C 

130=C 

190=  PROGRAM  RLCCHG ( INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT, TAPES) 

200=  DIMENSION  Y(20) ,X(20) ,GA(20) ,PE(20) ,XX(4 ,20) ,H(9) 

210=  DIMENSION  ZZ1 ( 20) ,ZZ2 (20) ,ZZ3( 20 ) , I FL( 7 ) 

220=  REAL  I,L,LA 

230=C  SECTION  1  INSERTIONS  ******************************** 

240=  R=0.445 

250=  L=0.353E-3 

260=  C=147.17E-6 

270=  PI=3. 1415927 

280=  EA=250.0 

290=  A=R/(2.0*L) 

300=  OMECA0=1.0/SQRT(L*C) 

310=  OMEGAA=SQRT(  OMEGAO**2-A**2) 

320=C  IPMT=0  HERE  IF  SUMMARY  CURVES  ARE  DESIRED 

330=C  NRUN- 

340=  NRUN-l 

350=  DO  6  IRUN=1 ,NRUN 

360=  IPNT=  0 

370=  WRITE(6 ,100) 

380=  I STOP-O 

390“  KEEP=1 

400=  IFL(1)=-1 

410-C  NIMT-  NUMBER  OF  INTEGRATORS 

420=  NINT=2 

430=  1  IFL(2)=— 1 

440=  2  CALL  RKFOUR(f'TINT,Y,X,GA,PE,XX,H,ZZl  ,ZZ2,ZZ3,TIME,PNTS,XEEP, 

IFL) 

450=  IF(IFL(4) )3 ,4 ,5 

460=  3  lriIT-1 

470“  TIME1=-1. 256739 

480=  NC=1 

490=  MCC= 1 

500-C  NCC=NUMBER  OF  POINTS  KEPT 

510-C  PINTIM-FIMISH  TIME 

520=  FINTIM*PI/OMEGAA 

530-C  PT1AX= MAXIMUM  NUMBER  OF  POINTS 

540=  PTMAX=800 
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550—C  SECTION  2  INSERTIONS  ************************************ 

560=  Y(l)=-5.4 

570=  H(l)=1.0E-07 

580=  GA(l)=1.0/C 

590=  GA(2)=1.0/L 

600=  WRITE(6 , 101 ) 

610=  4  IF(XIME-TIM£1)7,8,7 

620=  7  TIMEl=TIME 

630=C  SECTION  3  INSERTIONS  **************************** 

640=  V=F,A*SIt!(OtlF.GAA*TI!lE) 

650=  8  CONTI MO E 

660=C  SECTION  4  INSERTIONS  ******************************* 

670=  VR=Y(2)*R 

680=  VL=V-VR~Y(1) 

690=  X(2)=VL 

700=  X( 1 )=Y(2) 

710=  I=Y(2) 

720=  VC=Y( 1 ) 

730=  GO  TO  2 

740=  5  IPR=0 

750=  NC=NC-l 

760=  IF(NC)9 ,9,10 

770=  9  NC=NCC 

780=  IPR=1 

790=  10  IF(TIME.CE.FINTIM.0R.PNTS.GF,.PTMAX)IST0P=1 
800=  IF((IPR+ISTO?).EQ.O)GO  TO  11 

810=C  SECTION  5  INSERTIONS  ******************************* 

820=  URITE(8 , 102) TIME, I, VC, VL 

839=C  URITE(6,102)TIME,I,VC,VL 

840=  11  INIT=0 

850=  IF(IST0P)1 ,1,6 

860=  6  CONTINUE 

870=  WRITE (6, 100) 

880=  100  FORMAT( 1H1 , IX) 

890=  101  FORMAT(7X,4HTIME, 13X,2HIC, 13X,2HVC, 13X,2HVL) 

900=  102  FORMAT(2X,E12 .5 , 3X,F,12 .5 ,3X,E12 .5  ,3X,E12 .5 ) 

910=  END  FILE  8 

920=  STOP 

930=  END 


PROGRAM  PJLCDIS 


IOO=C  THIS  PROGRAM  USES  THE  RKFOUR  SUBROUTINE  TO  INTEGRATE  THE 
110=C  DIFFERENTIAL  EQUATIONS  IN  AN  P.LC  DISCHARGE  LOOP.  THE  SOURCE 
120=C  VOLTAGE  IS  SET  TO  ZERO  AND  THE  CAPACITOR  IS  CHARGED  TO  AN 
130=C  INITIAL  VALUE.  THE  INITIAL  CURRENT  IS  SET  TO  ZERO. 

140=C  THE  OUTPUT  IS  READ  TO  TAPES  WHICH  MAY  THEN  BE  EDITED  FOR 

150=C  GRAPHING  OR  VIEWING.  THE  TIME  AT  WHICH  INTEGRATION  CEASES  IS 

160=C  DEPENDENT  UPON  THE  DAMPED  RESONANT  FREQUENCY  OF  THE  CIRCUIT. 

170=C 

180-C 

190=C 

200=  PROGRAM  RLCDIS( INPUT ,0'JTPUT, TAPE5=INPUT,TAPF.6=OUTPUT, TAPES ) 

210=  DIMENSION  Y( 20) ,X(20) ,CA(20) , PE(20 ) ,XX(4 , 20) , H( 9) 

220=  DIMENSION  ZZ1 ( 20 ) ,ZZ2(20) ,ZZ3( 20) , IFL( 7 ) 

230=  REAL  I.L 

240=C  SECTION  1  INSERTIONS  ******************************** 

250=  R-II.O 

260=  L=29.0E-6 

270=  C=0.43E-6 

230=  PI=3. 1415927 

290=  A=R/(2.0*L) 

300=  OMECAO=1 .0/SQRT(L*C) 

310=  OMEGAA=SQRT( OMEGAO**2-A**2) 

320=C  IPNT=0  HERE  IF  SUMMARY  CURVES  ARE  DESIRED 

330=C  NRUN= 

340=  NRUN= 1 

353=  do  6  iru::=i  ,  "j:: 

36'.:=  ip:;t=  o 

370=  ITF.(6 , 100 ) 

3-'.0=  IST'.v>=0 

390=  XF.KP-1 

400=  IFL( 1 )=-l 

4io=c  mi::t=  number  of  integrators 

420=  NINT=2 

430=  1  IFL(2)=-1 

440=  2  CALL  RKFOUR(NINT, Y,X,GA,PE, XX,  11, ZZ1 .ZZ2.ZZ3 .TIME, PUTS .KEEP, 
450=  *IFL) 

460=  IF(IFL(4))3,4,5 

470=  3  INIT= 1 

480=  TIMEl=-l. 256739 

490=  NC= 1 

500=  NCC=1 

5 10=C  NCC=NUMBER  OF  POINTS  KEPT 

520=C  FIMTIM=FINISH  TIME 

530=  FIMTIM=Pl/OMEGAA 

540=C  PTMAX= MAXIMUM  HUMBER  OF  POINTS 

550=  PTMAX=800 


SECTION  2  INSERTIONS  ******************************* 
570=  Y(l) =-6000.0 

580=  H( 1 )=1 .0E-07 

590=  GA(1)=1.0/C 

600=  GA(2)=1.0/L 

610=  WRITE(6 ,101) 

620=  4  IF(  TIME-TIMED  7, 8, 7 

630=  7  TIME1=TIME 

640=C  SECTION  3  INSERTIONS  **************************** 

650=  V=0.0 

660=  3  CONTINUE 

670=C  SECTION  4  INSERTIONS  ******************************* 
680=  VR=Y(2)*R 

690=  VL=V-VR-Y(1) 

700=  X(2)=VL 

710=  X( 1 ) =Y ( 2 ) 

720=  I=Y(2) 

730=  VC=Y( l ) 

740=  GO  TO  2 

750=  5  IPR=0 

760=  NC=NC-1 

770=  IF(NC)9 ,9 , 10 

780=  9  NC=HCC 

790=  IPR=1 

800=  10  IF(TIME.GE.FINTIM.OR.PMTS.GE.PTMAX)ISTOP=1 
810=  IF( ( IPR+ISTOP) .EQ.O)GO  TO  11 

320=C  SECTION  5  INSERTIONS  ******************************* 

830=  V/RITE(8 , 102)  TIME  ,I,VC,VL 

840=C  URITE( 6 , 102 )TIME ,  I , VC , VL 

850=C  CALL  FPLOT(800,IPNT,AR,LR,0,1 ,2, TIME, V) 

860=C  CALL  FP LOT (800 , IPNT ,AR,LR, ISTOP ,2 , 2 ,TIME , VC) 

870=  11  INIT=0 

880=  IF( ISTOP ) 1,1,6 

890=  6  CONTINUE 

900=  V;RITE(6,100) 

910=  100  FORMAT(  1111 ,  IX) 

920=  101  FORMAT(7X,4HTIME,13X,2HIC,  13X,2I!VC,13X,2I1VL) 

930=  102  FORMAT(2X,E12.5,3X,E12 .5 ,3X,E12.5,3X,E12.5) 

940=  END  FILE  8 

950=  STOP 

960=  END 


Appendix  D 


Equipment  List 

Alternator  -  3  phase,  A  wire,  24  pole 
Model  LPL-72-2450 
S/N  711941-0007 

18.8  KVA,  15  KW  @  0.8  pf  lagging 

120/208  volts,  52  amps 

240/416  volts,  26  amps 

400  Hertz,  2000  RPM 

=  0.0594  pu 

X  ,  *  0.456  pu 
ad 

X  =  0.265  pu 

aq 

High  Leakage  Transformers  (3) 

Electro  Type  E  18230;  S/N  1>  2,  and  3 
19.98  KVA,  400  Hertz,  Single  phase 
Terminals  1-2:  240  volts,  84  amps 

Terminals  3-4:  4440  NL,  4.5  amps 

DC  power  Supply  (Exciter  M  citation) 
HP  Model  62928 


Instrumentation 

,834  Storage  Oscilloscope  -  T.ktroul*  «lth  »V  probes 
Ammeters  -  Weston  Model  433 

Wattmeters  -  Westinghouse  S/H  32982/32983/32984 
Digital  Multimeter  (True  8KS)  -  HP  Mraiel  3466A,  S,»  2488S 
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Thyratrons  -  JAN  8613,  5  MW  Peak  Power 
20  KV  Peak  Anode  Voltage 
500  Amps  Peak  Forward  Current 


Charging  Diodes  -  R«tifi.r  Stork  Built  Bp  of  MK  of  resistors 
capacitors ,  and  1  KV  PRV  diodes 
30  KV  Peak  Reverse  Voltage 
20  Amps  Nominal  Current 


0.05  yf 
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Appendix  E 


Leakage  Transformer  Short-Circuit  Test 
(C  Phase  Transformer  S/N  3) 

Ip  =81.6  amps  RMS 

Ig  =4.5  amps  RMS 

Vp  =  142  volts  RMS 

P  =  950  watts 

R  =  P/Ip2  =  0.1427  ohms 

Z  =  Vp/Ip  =  1.74  ohms 

X  =  SQRT  (Z2  -  R2)  =  1.734  ohms 

L  =  X/[27t(400)  ]  =  0.69  mh 

Zfi  =  240/84  =  2.857 

Rpu  =  R/ZB  =  0.04995  pu 

XpU  =  X/ZB  =  0.6069  pu 
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